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Abstract: We introduce ultradiscrete tau functions associated with rigged 
configurations for . They satisfy an ultradiscrete version of the Hirota 
bilinear equation and play a role analogous to a corner transfer matrix for the 
box-ball system. As an application, we establish a piecewise linear formula for 
the Kerov-Kirillov-Reshetikhin bijection in the combinatorial Bethe ansatz. 
They also lead to general 7V-soliton solutions of the box-ball system. 

1. Introduction 

The Bethe ansatz and the corner transfer matrix are methods of primary importance 
in analysing solvable lattice models [1]. The Bethe ansatz produces eigenvectors of row 
transfer matrices from solutions of the Bethe equation [2]. The corner transfer matrix 
method determines the one-point function from the one-dimensional sums [1]. See [3, 4, 5] 
and [6, 7] for some typical applications. Interestingly, both of these approaches are known 
to admit combinatorial versions, which have brought fruitful insights and applications 
into representation theory as well [8] . 

The combinatorial Bethe ansatz was initiated by Kerov, Kirillov and Reshetikhin 
(KKR) [9, 10]. They invented the object called rigged configuration, which serves as a 
combinatorial substitute for the solutions of the Bethe equation. By the KKR bijec- 
tion, they are in one-to-one correspondence with the Littlewood-Richardson tableaux, or 
equivalently, highest paths which are the combinatorial analogues of the Bethe eigenvec- 
tors. As for the corner transfer matrix method, a decisive progress came with the advent 
of the crystal base theory [11, 12], where the one-dimensional sums are formulated as 
generating functions of the energy of affine crystals over paths. 

Guided by a number of relevant results [13, 14, 15, 16, 17, 18], these streams have 
merged into the so-called X = M conjecture [19, 20] for general affine Lie algebra. Here 
X is the one-dimensional sum in the corner transfer matrix method. For type A„ , 
it coincides essentially with the Kostka-Foulkes polynomial [21] for the case treated in 
[9, 10]. On the other hand, M is the fermionic formula (2.10) in the Bethe ansatz, 
which is a generating function of the charge function c(/x, r) (2.9). By now, the X = M 
conjecture has been studied extensively and solved in several cases [22, 23, 24, 25]. 

During these developments, it was realized that not only the Bethe ansatz or the corner 
transfer matrix, but also the solvable lattice models themselves admit decent combinato- 
rial versions. In fact, vertex models with the quantum group symmetry U q (A^') turned 
out to be the soliton cellular automata at q = [26, 27] that had been known as the 
box-ball systems [28, 29]. Row transfer matrices in the former tend to commuting time 
evolutions in the latter. The finding has led to a systematic generalization of such au- 
tomata [30, 31, 32], which possess fascinating features as ultradiscrete integrable systems 
[33]. (See the explanation under (5.10) for the ultradiscretization.) Thus it is a natural 
endeavor to study these automata by the combinatorial versions of the Bethe ansatz and 
the corner transfer matrix. 

As for the Bethe ansatz, this has been done in [34, 35], which yielded the inverse 
scattering formalism of the box-ball systems. It turned out that rigged configurations are 
action-angle variables, which provide the conserved quantities or linearize the commuting 
time evolutions. The KKR bijection is the direct/inverse scattering (Gel'fand-Levitan) 
map. In particular, the mysterious combinatorial algorithm in the bijection is identified 
with a crystal theoretical vertex operator. 

Then what about the corner transfer matrix? And this is the issue that we are going 
to address in this paper. From a naive point of view, one is tempted to regard the 
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number of balls in a quadrant of the two-dimensional time evolution pattern of the box- 
ball system as its candidate. We introduce such a quantity Pi{p) (4.1) for a path p. On 
the other hand, the combinatorial analogue of the corner transfer matrix in the crystal 
base theory is the energy of affine crystals [12, 17], which is denoted by £i{p) in (4.12). 
Our Proposition 4.6 asserts Pi(p) = £i(p) indeed. One of the main results in this paper 
is Theorem 6.12, which states n(p) — Pi(p) — £i(p)- Here n(p) is the piecewise linear 
function on the rigged configuration (//°), (p^ 1 \r^), . . . , r^)) for p: 

Ti(p) = max{-c(i/,s) - |^ w |}, 

c ^ s ) = \Y. min (^ (a) . vih) ) - min ( A > ) + E i s(a) i . 

a,b a 

where {C a b)i< a ,b<n is the Cartan matrix of A n . c(v, s) is the charge function appearing in 
the fcrmionic formula, and the max extends over all the subsets (^ a \s (a) ) C (p {a \r^) 
of the rigged configuration. See (2.19), (2.20), (2.24) and Section 2.1 for a precise account. 
In short, Tj is an ultradiscretization of a single summand in the fermionic formula with 
respect to the subsets of the rigged configuration. 

An origin of this curious quantity goes back to Sato's theory of soliton equations [36] . 
In fact, tj arises as an ultradiscretization of the well known tau function for the KP 
hierarchy [37] under a special choice of parameters adapted to the rigged configuration. 
Using this fact, we show that n satisfies an ultradiscrete version of the Hirota bilinear 
equation, which actually serves as a characterization of n up to a boundary condition. We 
call Tj the ultradiscrete tau function. It serves as an analogue of a corner transfer matrix 
in the box-ball system and bilinearize the dynamics. These features are summarized in 
the following table. 





Bethe ansatz 


Corner transfer matrix 


main combinatorial object 


rigged configuration 


energy in affine crystal 


role in box-ball system 


action-angle variable 


tau function 


description of dynamics 


linear 


bilinear 



As the main consequences of Theorem 6.12, we derive a piecewise linear formula for 
the KKR bijection (Theorem 2.1), the solution of the initial value problem (Theorem 7.6) 
and the general N-soliton solution (7.21), (7.37), (7.42) for the box-ball system. Note 
that the quantities pi — £{ arise from the corner transfer matrix and crystals, whereas Tj 
is an explicit formula originating in the Bethe ansatz. Therefore our Theorem 6.12, i.e., 
Pi = £i = Ti provides another connection of the two methods analogous to the X = M 
conjecture. 

The layout of the paper is as follows. In Section 2, is introduced in (2.18)-(2.20) as a 
piecewise linear function on rigged configurations. It is actually a member of the family 
(2.22) which obeys the recursion relation (2.23). It reflects the nested structure 
sl n+ i D sl n D • • • D sl 2 , which will be utilized extensively. The piecewise linear formula 
for the KKR bijection is stated in Theorem 2.1. 

In Section 3, we give the definition and the basic properties of the box-ball system. 

In Section 4, we introduce pi and £j. pi in (4.1) is the number of balls in the SW 
quadrant in the time evolution pattern of the box-ball system. £j is defined by (4.12) and 
(4.11), which is a sum of local energy function in the affine crystal. They are analogues of 
the corner transfer matrix [1] in complementary viewpoints; pi originates in the box-ball 
system and £i in the crystal base theory. They are identified in Proposition 4.6. 
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The piecewise linear formula for the KKR bijection (Theorem 2.1) is a consequence of 
the further identification r, = pi = £i in Theorem 6.12. Sections 5 and 6 are devoted to a 
proof of this fact. In Section 5, Ti is shown to emerge as an ultradiscretization of the tau 
functions of the KP hierarchy (Lemma 5.3) and satisfy the Hirota type bilinear equation 
(Proposition 5.1). The key to these results is the special choice of the parameters (5.5)- 
(5.9). It assures the positivity, which is vital in the ultradiscretization (Lemma 5.2). The 
content of this section is a refinement of the earlier analysis [26]. 

In Section 6, Tj = pi for An is proved on the asymptotic states by induction on 
the rank n (Proposition 6.1 and its reduction in Proposition 6.4). From the assumption 
Ti = Pi = £i for ^4^— l ' tnc scattering data is expressed in terms of tau functions (Lemma 
6.6). Then we take advantage of the vertex operator formulation of the KKR bijection 
[34, 35] to make the induction proceed. Combined with the results in Section 5, the 
agreement on the asymptotic states is enough to establish the claim n = p% everywhere. 

In Section 7, Theorem 2.1 and Theorem 6.12 are generalized to arbitrary (non-highest) 
states. As an application, we present the solution of the initial value problem of the 
box-ball system in Theorem 7.6. Our tau functions are parametrized by the conserved 
quantities that specify solitons. We rewrite them in several forms in (7.21), (7.37) and 
(7.42). They yield general TV-soliton solutions of the box-ball system. Among others, 
our ultradiscrete tau functions are most elegantly presented in (7.42) in terms of affinc 
crystals in the "principal picture" . 

Appendix A summarizes the rudiments of the crystal base theory. Appendix B il- 
lustrates the graphical rule [17] for obtaining the combinatorial R, the winding and the 
non-winding numbers relevant to the energy function. Appendix C recalls the combinato- 
rial algorithm for the KKR bijection. Appendix D is the crystal theoretical reformulation 
of the KKR map due to [34, 35] . Appendix E is an exposition of the inverse scattering 
formalism of the box-ball system which supplements Section 3. 

2. Ultradiscrete tau function 

2.1. Preliminary. We summarize the basic notation used throughout the paper. For a 
multiset A = (Ai, . . . , Afc), we use the symbols 

(2.1) |A| = Ai H h A fc , £(X) = k, 

(2.2) A[jv] = (Ai, . . . , Ajv), (0<N<k), 

where A[o] = 0. Given two multisets A = (Ai, . . . , Afc) and p = (pi, . . . , p m ), we use the 
notation: 

k m 

(2.3) min(A, p) = min(A^, pj), 

i=l j=l 

(2.4) AC/i {Ai, . . . , Afc} C {p 1: p m }, 

where C accounts the multiplicity as well. For example, 0, (1, 1), (1, 3, 1) C (1, 2, 1, 3) but 
(2,2) % (1,2,1,3). 

2.2. Rigged configurations. Consider the data of the form 

(2.5) ( M (°',(^,r (1, ),..,(^ ) ,r("))), 

where p^ = [p[ a) , . . . , p^) G (Z>i)'« and r< a ) = (rj a) , . . . , ) G (Z> )'» for some 
l a > 0. Apart from p(°\ each (p^ a \A a ^) is to be understood as a multiset of the pairs 
(p^i\ r[ a ^ ),..., (Pi^ , r^) whose ordering does not matter. The data (2.5) is called a 
rigged configuration for if 

(2.6) < 4 a) < p {a l } for any pair (^ a) ,r t (a) ). 
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Here p^ is called the vacancy number and defined by 

(2.7) = E { f- 1] - 2E { f ] + E^ a+1) (1 < a < n), 

(2.8) E\ a) = J2 min(j, p[ a) ) (0 < a < n), 25< n+1 > = 0. 

fc=i 

The array (fi^ , . . . , ) is called a configuration and the nonnegative integers r are 
called rigging. Note that p^ and E^ depend only on the configuration. In particular 
E^ = It is customary to arrange as fi^ = (fi^ > ■ ■ ■ > p,^) and regard 

the rigged configuration as an n-tuple of Young diagrams pS x \ . . . , p^ where the row 
of length jif" 1 is assigned with the rigging subject to the condition (2.6). In this 
convention, we identify all the diagrams obtained by reordering the rows of equal length 
with different rigging. In what follows we do not assume p!^ > ■ ■ ■ > p^ unless explicitly 
mentioned. 

For a multiset with positive components A, let RC(A) denote the set of rigged config- 
urations (2.5) with p^ = A. Set 

(2.9) c(p, r) = \ ]T C ab mm(p^ , p^) - mm(p^ , M «) + ]T |r<°> \ , 

a,b a 

where (C a b)i< a ,b<n is the Cartan matrix of A n . The fermionic formula [10, 9] is obtained 
as the generating function: 

(2.10) M(A) = ^2 q c ^ p ' r \ 

where the sum extends over all the rigged configurations (A, (p^ , r^), . . . , (p^ n \ r^™')) G 
RC(A) with prescribed values for l// 1 ^, . . . , \p^\- The sum (2.10) is arranged as M(A) = 

sum over the rigging r under the condition (2.6) yields 
a product of g-binomial coefficients as is well known. 

2.3. Crystals. We recapitulate basic facts on the crystal B\. For a general 
background see Appendix A. The B\ is the crystal base of the Z-fold symmetric tensor 
representation. As the set it is given by 

(2.11) B l = {x={x 1 ,..., x n+1 ) e (Z> )" +1 | &! + ■■■ + x n+1 = I}. 

The Kashiwara operators act as e^(x) = x' , fi(x) = x" with x'j = xj + 5ij — Sij+i and 
x" = Xj — Si j + Sij+i. Here indices are in Z n+ i and x' and x" are to be understood 
as unless they belong to (Z> )™ +1 . The combinatorial R : Aff(Bi) ® AS(B m ) — > 
AS(B m ) ® AS(Bi) has the form R : x[d]®y\e] i-> y[e - H(x®y)] ®x[d+ H(x ® y)}, which 
are described by the piecewise linear formula [38, 26]: 

(2.12) Xi = Xi + Qi(x (g) y) - Qi-i(x ® y), & = yi + Qi-i(x <S> y) - Qi(x <g> y), 

fc-1 n+l 

(2.13) Qi(x (g> y) = min{y^ x i+j + ^ Ui+j 1 < k < n + 1}; 

j=l j=k+l 

(2.14) H(x ® y) — min(Z,m) - Q (x ® y). 

The energy function H here is normalized so that < H < min(Z, m) and coincides with 
the "winding number" [17]. In general min(/,m) — Q t is the i th winding number that 
counts the lines crossing Xi and Xi + \ (Appendix B). 

The element x = (x\, . . . , x n+ \) is also denoted by a row shape semistandard tableau 
of length I containing the letter i Xi times and x[d] G Aff(£>;) by the tableau with index 
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d. For example in A^\ the following stand for the same relation under R: 
(1,2, 0, 1)[5] ® (1, 0, 1, 0)[9] ~ (0, 1, 0, 1)[8] ® (2, 1, 1, 0)[6], 

(2.15) 



1224 


5® 


13 


q — 


24 




1123 



To save the space we use the notation: 

(2.16) a 1 = \a---a\ E B[, u t = l' = |T- ■ ■ 1 | g B t . 
Setting 

B} a+1 = {{x u . . .,x n+1 ) e Bi | xx = ■ ■ ■ = x a = 0} (0 < a < n), 

we have 

(2.17) B t = B} 1 d Bp D ■ ■ ■ D B^ n+1 = {(n + l) 1 } 

as sets. We will need to consider the crystals not only for but also for the nested fam- 



ily Aq , A\ , . . . , A n _ 1 . In such a circumstance we realize the crystal Bi for A n _ a (0 < 



a < n) on the set Bf a+1 with the Kashiwara operators e^/i (a < i < n). In this 



convention the highest element with respect to A n - a is (a + 1)' e Bf a+1 . 
Let 

P+(A) = {p e B Xl ® • • • ® B At | = 0, 1 < i < n}. 

be the set of highest elements (paths) with respect to A n . The bijection [9, 10] between 
RC(A) and the Littlewood-Richardson tableaux is translated to the one between RC(A) 
and V+(X). We call the resulting map the KKR bijection. See Appendix C for an 
exposition of the algorithm and Appendix D for the recent reformulation as the crystal 
theoretical vertex operator [34, 35]. In particular, there is a nested structure with respect 
to the rank in the sense that if (n^ \ (^ 1 \r^), . . . , (^ n \r^)) is a rigged configura- 
tion for A [ n \ so is (fi (a \ (fi (a+1 \r (a +V), (//("), r< n ))) for A^ ] _ a . Moreover, the KKR 
bijection sends the latter to a highest path in B-f^f 1 ® • • • <g> B-^t 1 . 

2.4. Piecewise linear formula for KKR bijection. We use the notation defined in 
Section 2.1. Given a rigg ed configuration (^°\ r^), . . . , (/jW,rW)), we introduce 
the ultradiscrete tau functions t (A), ti(A), . . . , t„ + i(A) for A C //°) as follows: 

(2.18) t (A)=t„ +1 (A)-|A|, 

(2.19) T d {\) =max{-c(^,s)- \is (d) \} (1 < d < n + 1, \v {n+1) \ = 0), 

vC.fi 

(2.20) -c(v, s) = min(A, v {1) ) + m.m{v {1 \ v {2) ) + ■■■ + min(^ ( ™- 1) , v {n) ) 

- mm(v (1 \v {1) ) - min(i/ 2) , v {2) ) min(*/ ,l) , i/ n >) 

-| S W| | S (")|. 

In (2.19), max is taken over v — (v^ 1 ' , . . . , i/( n >), where the components are independently 
chosen under the condition z/ 1 ) C /i' 1 ), . . . , i/™) C //("). The array s = (s^ , . . . , s^) 
denotes the set of the riggings C , . . . , C that are paired with the chosen 
. . . , i» as {(^ (a) , s[ a) )} C {(^ (a) , ) )}. The quantity c(u, s) in (2.20) is obtained 
from c([i,r) (2.9) by replacing (n,r) = (//°\ . . . , (//"), r( n ))) with (^, s) = 

(A, (^' 1 - ) , s^ 1 )), . . . , (z^(™\ s^™^)). Apart from C ^' a \ there is no further constraint on 
|^W|,...,|z»| and it is not required that the data (A, (i/ 1 ), s^ 1 )), . . . , (v {n \ sW)) to be 
a rigged configuration for A„ . Since the max (2.19) includes the trivial case Mv^ = 
0, the quantities n(A), . . . ,r n +i(A) are nonnegative integers. Note that r„ + i(^°)) = 
max{— c(^, s)} in (2.19) may be viewed as an ultradiscretization of the single summand 
qc(p,r) - n f erm j on i c formula (2.10) with respect to the subsets (^, s) C (/x, r). See also 
(5.11). 



Theorem 2.1. Let the image of the rigged configuration (^°\ r^), . . . , {^ n \ r^)) 
under the KKR bijection be the highest path p\ ® • • • <g> Pl G V+{p^). Then p k = 
(xi, . . . , x n+ i) G B (o) is expressed as 

(2.21) X d = Tk,d - Tk-l,d - Tk,d-1 +Tk-l,d-l, 

where r M = T d ((fi { °\ . ..,(/° ) )). 

Note that (2.18) ensures x\ + ■ ■ ■ + x n+ \ = n k °\ 

Due to the nested structure of the KKR bijection with respect to the rank [34], The- 
orem 2.1 is also stated as a family of relations corresponding to sl n+ i D sl n D ■ ■ ■ D sfo- 
To do so, wc introduce the family of ultradiscrete tau functions {T d a \\) < a < 
n-1, a<d<n + l, A C ^} by T^ a) (A) = ^(A) - |A| and 

T { d a) (X) = max{min(A, !> +1 >) + min(i> +1 \ i/ a+2 )) + ■■■ + min^™- 1 ), i») 

(2.22) - min(i/° +1 >, - min(z> +2 \ i/ a+2 >) min(i/W, i/W) 

- | S ( a+1 )| - | S ( a+2 )| - |^|} (a + 1 < d < n + 1), 

where = as before. The max is taken over the independent choices j/ a+1 ) C 

^ a+1 ), . . . , i» C //"). The subsets of the riggings s< a+1 ) C r ( a+1 \ . . ., C r (") are 
those paired with the chosen before. The previously introduced tau 

function T d (X) (2.19) is equal to r^(A). Now Theorem 2.1 is rephrased as 

Theorem 2.2. Given a rigged configuration (fJ,^°\ (^ 1 \r^), . . . , r^ n ^)) and < 

a < n-1, let the image o/(/^W,( ^ l( a+ V ( '' +1) ),■■• I (/ 1, " ) .'' ( " , )) <fte & *" 

jection be the A n - a highest path p\ ® ■ ■ ■ ® pi £ B-.^f 1 eg) • • • B-.1t 1 . Then p k = 

i'\ I'),,' 

(x a+ i,x a+2 , ■ ■ ■ , x n+ \) is expressed as 

(a) (a) (a) (a) 

x d — T k,d ~ r fc-l,d - T fe,d-1 + r fe-l,d-l' 

«^rW=rW(( M i o) ,...,M£ ,) ))- 

Again, x a+ i + • • • + x n +\ = is evident by the construction. For a proof of Theorem 
2.1, see Section 4.4. 

The tau functions (2.22) are the solution of the recursion relation with respect to the 
rank: 

(2.23) rW(A) = max {min(A, v) - min(i/, v) \s\ + t<° +1 >(i/)} 

vC.fj,\ a + 1 ) 

for < a < n — 1, a+l<d<n+l with the convention ri a '(A) = t^+^A) — |A| and 

the initial condition t^™\(A) = 0, r^(A) = — |A|. The rigging s is the subset of r( a+1 ) 
paired with the chosen v. 

Lemma 2.3. T^ a '(0) = for any < a < n — 1 and a+l<d<n+l. 

Proof. It suffices to prove a — case. When A = 0, (2.19) becomes T d {%) = — min^c/i{c(f)+ 
|sW| + • • • + |s(")| + |s( d )|}, where c(v) is given by (see (2.9)) 



5 H = |E c ^ min ^ (a) ^ (b) ) = jE C ^E min ( l 'J> 



a, b a, 6 



where is the number of k such that i>j^ = i. This is a positive definite quadratic 
form whose minimum is at Vmj"' = 0. The other part \s^\ + • • • + \s^ \ + \s^\ 
appearing in 7^(0) also attains the minimum simultaneously at this point. □ 
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Let the image of the rigged configuration (^°\ (^ 1 \r^), . . . , {^ n \ r^)) under the 
KKR bijection be the highest path pi ® ■ ■ ■ ® Pl G P+(//°)) C B <o) <8> ■ ■ ■ <E> B ( >. In 

'Mi Mi 

what follows we will also write 

(2.24) T i (A)=T M = T i (pi®---®p fe ) for A = /ig = (a4°\ • • • , ) (1 < k < L). 

Concerning the notation Tj(pi (8) ••• <S> Pk), a remark is in order. Any highest path 
Pi <8> • • • <S> Pfe can be extended to a longer one pi <8> ■ ■ ■ ® Pfe ® Pfc+i ® • • • ® PL in 
which pfc+i ® ■ ■ ■ ® pl is not unique. Suppose that (//°), (/t^ 1 ), r^), . . . , (//™\ r^)) 
and (/i (^i W,r (- 1 )), . . . , (/i ("),r < -™' ) )) are two rigged configurations corresponding to 
such extensions of pi ® • • • <g> pfe, and let Tj(pi ® • • • <8> Pfe) and r t '(pi ® ■ • ■ ® Pk) be the 
associated tau functions in the sense of (2.24). Then r,(pi ® • • • ®Pfe) = r/(pi <g> • • • £g>Pk) 
will be guaranteed by Theorem 4.9. Note however that they are different as the piecewise 
linear expressions as in (2.18)-(2.20). By the reason, we will always mention the rigged 
configurations relevant to pi ® • • • ®pfc. 

Example 2.4. Consider the highest path p = 11112221322433 e Bf L of length L = 14, 
where we have omitted the symbol <g>. The corresponding rigged configuration is depicted 
in Example C.2. Thus we set 

^(l 14 ), M W = (4,3,2), M (2) = (3,l), ^ (3) = (1), 

r* 1 ) = (0,2,3), r< 2 ) = (l,0), = (0). 

The associated tau function Tk,i takes the following values. 



k 
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2 
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10 


11 


12 


13 


14 


Tk,l 
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Tk,2 
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11 


13 
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1~k,3 
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12 
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18 


7fc,4 
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6 


8 


10 


13 


16 


19 



The choices of the subsets v = (i/ 1 ) , v^ 2 \ v^ ) that attain these values for Tk,4 — 
maxj,c P {' • • } in (2.19) are as follows. 



k 


V 


1,2,3 


A 


4 


A, B 


5,6,7 


B 


8 


B,C 


9,10 


C 


11 


C,D 


12,13,14 


D 



Here A,B,C,D C fi = (/Z 1 ), // 2 ) , fj,^) are given by 
A = (0,0,0), 

B = ((4),0,0), ((4),(1),0), ((4),(1),(1)), 

C = ((4,2),(3),0), ((4,2),(1),0), ((4,2),(3),(1)), ((4, 2), (1), (1)), ((4, 2), (3, 1), (1)), 
C = ((4,3,2),(3,1),(1)) = ( A1 ( 1 ), M ( 2 ), M ( 3 )). 

The case k = enforces the choice \/v^ = in agreement with Lemma 2.3. In the other 
extreme case k — L, the full choice v = fi is the consequence of the general result in 
Remark 6.14. In general the maximum attaining v for r»(A) = max„c/j{' ■ ■ } gradually 
grows with A. The above p will be investigated further in Examples E.l and E.4. 
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3. Box-ball system 

3.1. Conventional formulation. Consider the tensor product B\ ± ® B\ 2 <g> • • • ® Saz,- 
Its elements are called states. We regard each component (xi, . . . , x n +\) G Bi as a 
capacity / box containing ajj balls with color i for 2 < i < n + 1 . On the other hand the 
letter 1 is to be interpreted as a vacancy. Thus x\ represents the empty space in the box. 
A state represents an array of boxes with capacity Ai, . . . , Xl containing balls of colors 
2,3,...,n+l. 

We define the time evolution Ti (p) = p' x ® • • • <g) p' L of a state p = p\ ® • • • ® pl by 

(3.1) ui [0] <g> Pl [0] <g> • • • ® PL [0] ~ P [ [-d x ] ® • • • ® p' L [~d L ] ® vi [di + • • • + d L ] 

under the isomorphism Aff(B ; )®(Aff(B Al )(X)- • -®AfF(.B A J) ~ (Aff(B Al )<S>- • -<gAfF(.B Aj ,))<S> 
Aff(_B;). Here vi G £>; and dj are uniquely determined by (2.12)-(2.14). We set 

(3.2) Ei(p) = ei-\ he L , ej — min(A.,, I) - dj, 

which has the property Ei{p®Uk) = Ei(p) for any k and 

It is known [31, 27, 26] that 7] is weight preserving, the commutativity XJTfc = T^T; is 
valid and Ei(p) is a conserved quantity, i.e., Ei(Tk(p)) = Ei(p) for any k and I, provided 
that pj = u\ j for L' < j < L with sufficiently large L — L' . The proof of these facts 
is based on the Yang-Baxter equation of the combinatorial R (Proposition A.l) and the 
property: 

(3.3) vi = ui if pj = u\ j for L' < j < L with sufficiently large L — L'. 

Ti stabilizes for I ^> 1, which will be denoted by T^. 

Since each dj is the winding number (2.14), Ei(p) is the sum of the non- winding 
number ej. In particular for I = oo, ej is equal to the number of balls Xi + • • • + x n+ \ in 
the j th box pj = (xi, x 2 , ■ ■ ■ , x n+1 ) G B\ j . Therefore we find 

(3.4) E oa (p) = number of balls contained in p. 

In the terminology of solvable lattice models, Ei is the energy associated with a row 
transfer matrix. It should not be confused with another energy £i (4.12) relevant to the 
corner transfer matrix. Their relation is given in Proposition 4.8. The conserved quantity 
Ei will be evaluated explicitly for highest states in Proposition 6.15 and for general states 
in Proposition 7.7. 

Example 3.1. The time evolution of the top row p under X^, i.e., T 00 (p), T^ (p), 
T^(p) are listed downward. The frame of the semistandard tableaux and the symbol ® 
are omitted. 

11 122 2 1333 1 1 4 1 1 1 1 1 1 1 1 1 1 
11 111 1 1222 3334111111111 
11 111 1 1111 2223433111111 
11 111 1 1111 1112322433111 

The conserved quantities are given by Ei(p) = 3, E 2 (p) — 5 and Ei(p) = 7 for I > 3. 

The time evolution can be calculated by a simple prescription [26]. We introduce 
a map Li (2 < i < n + 1) by 

L t : Z> x Bi — > B[ x Z> 
( - 3 ' 5 ' ) {m,y)i — >(y',m'), 

where m' and y' = (y[ , . . . , y' n+ i) are determined from m and y = (yi, . . . , y n +i) by 

!yi + (yi-m)+ if .7 = 1, 
min(m, j/i) H j = i, 

yj otherwise, 
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where (m) + = max(m, 0). Lj may be viewed as the interaction of the box B t with the 
carrier that contains m balls of color i. The carrier drops as many balls as possible into 
the empty space y\ and picks away all the color i balls that were originally in the box. 
Using Li, we introduce the operators Ki (2 < % < n + 1) that sends a state to another as 
follows. 

Ki{p\ <g)f>2 ® • • • ) = Pi ® P2 ® • • ■ , 

Li{{mj,pj)) = {pj,m j+1 ) for j > 0, (ra = 0). 

The latter relation is applied successively for j = 0, 1, 2, . . ., determining all the p^-'s. In 
other words the operator Ki attaches an empty carrier to the left of the state and sends 
it to the right, by which the color i balls are moved to the right according to the local 
interaction rule Li. 

Proposition 3.2 ([26]). The time evolution admits the factorization: 

Too = K 2 K 3 ■ ■ ■ K n+ i. 

Example 3.3. For p in Example 3.1, K 4 (p), K 3 K 4 (p) and K 2 K 3 K 4 (p) = T^ip) are 
given. 

11 122 2 1333 1 1 4 1 1 1 1 1 1 1 1 1 1 
11 122 2 1333 1 1 1 4 1 1 1 1 1 1 1 1 1 
11 122 2 1111 3334111111111 
11 111 1 1222 3334111111111 

Remark 3.4. Suppose pj = u\ } for 1 < j < k in a state p = p\ ® ■ ■ ■ ® pl- Then 
Proposition 3.2 tells that in the state T 00 (p) = p[ ® • • • ® p' L , p'- = u\ j is valid for 
1 < J < fc + 1 . 

3.2. Bethe ansatz. Highest states in B\ x ® • • • ® B\ L are in one to one correspondence 
with rigged configurations {^°\ (fi^, ) , . . . , (fi^ n \ r^ )) with (j,^ — (Ai, . . . , Al) by 
the KKR bijection. Suppose L is sufficiently large. If a state p = p\ ® ■ ■ ■ ® pl is highest 
and pk — u\ k for k 3> 1, so is its time evolution Ti(p). Thus the box-ball system induces 
the time evolution on the associated rigged configurations. For such states, (2.8) and 
the vacancy number p^p are sufficiently large, and one can increase the color 1 rigging 
r- 1 ' without violating the condition (2.6). 

Proposition 3.5 ([34], Proposition 2.6). Let p = pi® ■ ■ ■ ®pl & V+{^^) be the image 
of the rigged configuration (fi^°\ (m^\ r^), . . . , (^ n \ r^™')) under the KKR bijection. 
Assume that vi = ui in (3.1) and set r/ 1 ^ = r^P + min(Z, fJ-P)- 

Then (fJ-(°\ {/J^, r (fiW , r^), . . . , (fJ^ n \ r^)) is a rigged configuration and corre- 
sponds to the highest state Ti(p) G V+{ji^). 

This is proved from the definition of the time evolution (3.1) and Lemma C.3. The 
time evolution TJ in this paper corresponds to the a — 1 case of T^ considered in [34] . In 
this sense the rigged configurations are the action-angle variables of the box-ball system 
which linearize the original nonlinear dynamics (3.1). Moreover it is clear that all the 
Ti(p) are the same if I > max^' 1 '. 

Example 3.6. The rigged configuration corresponding to T^(p) (t = 0, 1,2,3) in Ex- 
ample 3.1 (apart from n^). 



Gu«r (1) ) (^ ( V (2) ) (n 



(3) r (3) 



12 








3t 









12 








3t 







n 


13 




3 + t 





Do 
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The length of each row is pf' and the numbers on its right and left are the rigging r\ a ' 
and the vacancy number p^ a) , respectively. (Vacancy numbers are exhibited here for a 
check of (2.6).) 

The Bethe ansatz produces transfer matrix eigenvectors from solutions to Bethe equa- 
tions. The KKR bijection is its combinatorial version in the sense that the former is 
replaced by highest states and the latter by rigged configurations. Thus we see that the 
combinatorial Bethe ansatz provides a linearization scheme, or equivalently, an inverse 
scattering method of the box-ball system [34] . See Appendix E for a further exposition 
combined with the vertex operator formalism of the KKR bijection. 

4. Corner transfer matrix 

4.1. Number of balls in the SW quadrant. Let p = p\ (g>- • -®>pl be a state and write 
its time evolution as T^(pi ® • • -®Pl) = p\ <S> ■ ■ ■'S>P t L , with p* = (a;*- 1; £*- 2 , . . . ,a;*, n+1 ) G 
B\ j . We do not assume that p is highest. For < k < L and l<d<n+l,we define 
the function p k ,d{p) e z >o by (po,d(p) = °) 

k k 

(4.i) PkAp) = 1>?,2 + • • • + x% A ) + 1>$,2 + • • • + 4„ + l)- 

J=l t>l j=l 

Here the second term is finite due to Remark 3.4. In fact the double sum may well 
be replaced by X)t=i ^j— t+i onr y where the nonzero contributions are contained. This 
region is depicted as the SW quadrant of the time evolution pattern like Example 3.1. 




The first term in (4.1) is the number of balls of color 2,3, ... ,d contained in the top row, 
which is the truncation p\ ® ■ • ■ ® pk of the state p. The second term counts the balls of 
all colors 2, . . . , n + 1 within the hatched domain. By the definition, pk.n+i is the total 
number of balls within p\ ® • • • <g) pk and the SW quadrant beneath it. Thus 

(4-2) Pfe,l(P) = Pk,n+l{Too{p)) 

holds. Note that p ky d(p) is independent of pk+i,Pk+2, ■ ■ ■ ,Pl- In this regard, we will also 
use the notation 

(4.3) p d (pi <g> • • • ® p k ) = p k ,d(p)- 
From Remark 3.4 it follows that 

(4.4) p d (m ®p\®---®Pk)=Pd{P\®---®Pk) 
for any I. 

The above picture reminds us of Baxter's corner transfer matrix (CTM) in solvable 
lattice models [1]. In fact pk y d serves its ultradiscrete analogue adapted to the box-ball 
system as we will see below. 
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Example 4.1. For p in Example 3.1, Pk,d{p) takes the following values. 
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19 



4.2. Bilinearization of box-ball system. By the definition, the k th component p^ = 



{x\, . . . , x n+ i) G B\ k in a state p ~ p\ ® ■ ■ ■ ® pl is expressed as 

(4.5) x d = p k ,d ~ Pk-i,d ~ Pk,d-i + Pk-i,d-i (l<d<n+l), 
where pk,d — Pk,d{p) for 1 < d < n + 1 and the extra one Pk,o(p) is specified by 

(4.6) Pk,o(j>) = Pk,n+i(j>) - (Ai H h A fe ) 



so as to satisfy x± + ■ ■ ■ + x n+ i = Afe. The formula (4.5) may be viewed, in a certain 
sense, as an ultradiscrete analogue of the Baxter formula (eq. (13. 1.12) in [1]): (cri) = 
Tr(SABCD)/Tv(ABCD) for one point function in terms of CTMs. 
We use the notation 

(4-7) p k ,d = Pk,d(Too(p)). 

Thus (4.2) reads 

(4-8) Pk,l =Pfe,n+l- 

Proposition 4.2. For 2 < d < n + 1 the following relation holds: 

(4-9) Pk,d-i + Pk-i,d = max(7ife !d + pk-\,d-i, p k -i,d-i + Pk,d - Afe). 

A similar fact has been shown in [26] . 

Proof. In the time evolution = K2K3 ■ ■ ■ K n+ i (Proposition 3.2), let us calculate the 
effect of the operator Kd on the k th boxpfe = (x\, . . . , x n+1 ) G B\ k in Kd+\ ■ ■ ■ K n+1 {p). 
In the following, the fact that color d balls are touched only by K& is taken into account. 
Suppose that the carrier contains m and m! balls with color d just before and after the 
interaction Ld (3.5). In (3.6) we are to set 

k 

m! = ^2{pj,d - Pj-i,d - P 3 .d-i + pj-i,d-i) - (p->p) 
3=1 

= (Pk,d - Pk,d-l) - {Pk,d - Pfc,d-l); 

m = m'|fe^fe_i, 

Vd = Xd = pk,d - Pk-\,d - Pk,d-\ + Pk-\,d-l, 

where we have used (4.5). As for the empty space yi concerning Ld in (3.5), we show 
that it is given by 

(4-10) yi=^k+P k ,d-Pk-i,d-Pk,d + Pk-i,d (2<d<n + l) 

by induction on d in the decreasing order d = n + 1, n, . . . , 2. In so doing, the bilinear 
relation (4.9) will be established simultaneously. 

For d = n + 1, (4.10) coincides with x\ in (4.5) by (4.6) and (4.8), hence it is correct. 
Then the relation m! = yd + (m — yi)+ (3.6) leads to (4.9). The new empty space is 
determined from y[ = yd + (yi — m) + = m! + y\ — m and is equal to 

Afc + Pk,d-i - Pk-i,d-i - PkA-l + Pk-l,d-l- 
This coincides with (4.10) with d replaced by d — 1, making the induction proceed. □ 
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The relation (4.9) is an ultradiscrete analogue of the Hirota bilinear equation. In view 
of (4.8), it determines pk-i,i,pk-i,2, ■ ■ ■ ,Pk-i,n+i successively from {p fe _ M , p M , pk,d 
1 < d < n + 1}. Thus all the pfc,d(^oo(p)) are fixed uniquely from the data at sufficiently 
large t and k. Then the local states are specified by (4.5). In this sense the ultradiscrete 
CTM pd achieves a bilinearization of the dynamics of the box-ball system. 

4.3. Relation to energy function. Let p G B\ 1 ® ■ ■ ■ ® B\ L be any element which is 
not necessarily highest. For 1 < k < L, we introduce the sum: 

(4.11) £nPi®---®Pk)= ]T Q^(P3®P^ +1) ) (l<i<« + l), 

l<j<m<k 

where Qi is the i th non-winding number (2.13) with the convention Q n+ \ = Qq. The 
element defined by sending p m to the left by applying the combinatorial R 

successively as 

Pj ®p 3+ \ ® • • • ®p m -i ®p m ^Pj ®Pj+i ® ■ ■ ■ ®P^ ri " 1) ®p' m ~i — 

^ Pj ® P J+ i ® pii +2) ® ■ ■ ■ ® p'^ 

- Pj ® P ( A +1) ® P'j+i ®---®p' m -i- 
We understand that (4.11) is for k = 0, 1. Using £V we define the i th energy Ei by 

(4.12) £ l {p\®---®Pk) = £ y l (u oa ®p x ®---®p k ) (l<i<n + l), 

where Uoo actually means u\ with sufficiently large I. £, does not depend on such I. In 
fact, from the graphical rule in Appendix B, we find Qi (ui ® x) = X2 + + • • • + Xi if x = 
(xi, . . . , x n+ i) and x\ + ■ ■ ■ + x n+ i < I is satisfied. Thus writing p^p = (xj t i, . . . , Xj in +{), 

(4.12) is split into the boundary and the bulk parts as 

k 

(4.13) £i(pi ®---®p k ) = ^2(x h2 + ■■■ + Xj.i) + £?{ Pl ® ■ ■ -®pk). 

3=1 

In particular, one has £i{p\ ® ■ ■ ■ ®Pk) = '{pi ® ■ • -®Pk) ifpi®---®Pk is highest. We 
warn that the quantity usually called energy [17, 20] is — £ n +\ or — £^+i up to an additive 
constant. In what follows, whenever the notation is used, it should be understood as 
m with sufficiently large I and the relevant quantity is independent of such I. 
To the relation x ® y ~ y ® x with e = Qi(x ® y), we assign the diagram 

x y 

(4.14) /\ 

y x 

where the suppressed i is to be mentioned nearby if necessary. 

Let a((xi, X2, ■ ■ ■ , x n+ \)) — (x 2 , Xs, . . . , xi) be the Dynkin diagram automorphism 
acting on Bi decreasing the tableau letters cyclically by one. We extend it naturally to 
the tensor product by a(p\ ® ■ ■ ■ ®Pk) = &{pi) ® ■ • ■ ® a{Pk)- Since the combinatorial R 
commutes with a, the i th non-winding number has the properties similar to the i = 
case. In particular, under the Yang-Baxter relation 

the equalities a + b = e + f and b + c = d + e hold. In fact, suppose the figure corresponds 
to AS(B k ) ® Aft(Bi) ® AS(B m ) -» AS(B m ) ® AS(B t ) ® AS(B k ) for some k,l and m. 
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If i = for instance, the associated non- winding number Q is related to H via (2.14), 
therefore by setting a = min(/c, I) — a, b = min(fc, m) — and c = min(Z, m) — c, the left 
hand side represents the following relation under the combinatorial R: 

z[Ci] ® 2/[C2] ® z[Cs] l/'tCa - a] ® z'[Ci + a] ® «[Cs] 

- l/'[C2 - a] ® z'[C3 - 6] ® i»"[Ci + o + 6] 

~ z"[C 3 - b - c] ® y"[(2 - a + c] ® x"[£i +0 + 6]. 

Similarly, by setting d = min(Z,m) — d, e = min(fc,ra) — e and / = min(fc,Z) — /, the 
same element is transformed along the right hand side as 

x[Ci] ® 1/K2] ® ^ s[Ci] ® ^*[Cs - d] ® y*[C 2 + d] 

~ z**[C 3 - d - e] (g) x*[Ci + e] ® y*[( 2 + d] 

~ z** [Cs - d - e] ® y** [C2 + d - /] ® x** [Ci + e + /] . 

Since the Yang-Baxter relation is valid among the affine crystals, we obtain not only 
x" = x**,y" = y** and z" = z** but also b + c — d + e, a — c = / - J and a + b — e + f, 
which are equivalent to the two relations b + c = d + e and a + b = e + f. Note that 
a + b + c^e + f + d in general. 

Remark 4.3. The energy is invariant under any reordering of pi ® • • • ® pk by the 
combinatorial R. Namely, ® ■ • ■ ® Pfc) = £i(p[ ® ■ ■ ■ ® p' k ) and (pi ® • • • ® Pk) = 

£i {p'i ® • • • ® p'k) hold if f>i ® • • • ® Pk — Pi ® • • • ® Pk by the combinatorial R. For 
i = n+l this is essentially Proposition 3.9 in [20] and the general i case follows from the 
symmetry under a. 

Let us consider a particular diagram involving pi ® ■ ■ ■ ® Pk, which is illustrated for 
k = 2, 3, 4. The general case is similar. 

Pi P2 P3 Pi 




Incidentally, this kind of diagrams have been known as the half twist in the construction 
of link invariants [39] . 

Lemma 4.4. The energy £^ {pi ® • • • ® Pk) is the sum of the non-winding numbers Qi 
(as e in (4-H)) attached to all the vertices of the corresponding diagram for pi ® • • • ®Pk 
as above. 

Proof. For k = 2 it is obvious. We use the definition (4.12) and illustrate the induction 
step along the one from k — 3 to k — 4. 




By the induction assumption, the sum of three • is equal to Xa<j<m<3 QiiPj ®Pm +1 ')- 
Thus we are to verify ei+e 2 +e 3 = X)i<i<4 Qi(Pj®P4 +1 ^)- But the Yang-Baxter equation 



14 

shown above tells that ei+e 2 +e 3 = di+d 2 +d 3 , and furthermore, d% = Qi(p3®p^) 1 d 2 = 
Q l {p2®pT),d l = Q t (jp 1 ®pf ) ). □ 

Lemma 4.5. p,(pi <g> • ■ • <g> p/c) - Pi(T oc (pi <g> • • • ® p fe )) = ei + • • • + e fe , w/iere ej 's are 
the i th non-winding numbers specified by the following diagram: 




Proof. In terms of the notation in (4.1), the difference of pi is evaluated as 

k k 

5Z( X i,2 + ' ' ' + Xj t i) + + ' ' ' + x ],n+l)- 

3=1 3=1 

By using the graphical rule [17] explained in Appendix B, it is easy to show that the non- 
winding number Q t (2.13) is given by ej = (x° 2 H ^ x j,i) + ( x j.i+i H ^ x ),n+i)- ^ 

The main result in this subsection is the following, which identifies the ultradiscrete 
CTM pi (4.3) with the energy £ j that originates in the crystal theory. 

Proposition 4.6. pi(p\ • ■ • <E>Pk) = £i(pi <8> ■ ■ ■ <S>Pk) holds for any k and 1 < i < n+ 1. 

Proof. For T^p) with sufficiently large i, its leftmost fc components become ua x <8> ■ ■ ■ ® 
UA fc due to Remark 3.4. In this case the both pi and £j are obviously zero. Therefore it 
suffices to show 

p»(pi ® • • • ® Pfe) - p»(p'i ® • • • ® p^) = f»(pi ® • • • ® Pit) - £i(pi ® • • • ® p' k ), 
where p'\® ■ ■ ■ ®p' k = T oa (pi g) • • • ®pfc). We illustrate the proof for k = 3. From Lemma 
4.5, we are to show £i(pi ® P2 ® P3) = £j(pi ® p 2 ® P3) + e i + e 2 + 63. Recall that 
p'l ® • • • ® p' fc is determined by carrying Woo by the combinatorial R through pi ® • • • ® Pfc 
to the right as Moo ®Pi ® ■ ■ ■ ®Pfc — Pi ® • • • ® P^ ® (•)• Combining this with Lemma 4.4, 
one can depict the two sides as follows: 




£i (Pi ® P2 ® Ps) ^(pi p' 2 ® j/ 3 ) + ei + e 2 + e 3 



We are to check a + b + c = a' + b' + c' + e[ + e' 2 + e' 3 . From Remark 4.3, we may assume 
Ai > A2 > A3 without loss of generality. Then the above equality is a consequence of 
the separate ones e[ = 0, a = a' + e' 2 and b + c = b' + c' + e' 3 . To see them, note that 
Moo ® b ~ u m ® (•) for any b E B m under the combinatorial R. Moreover Qi(u m ® Uj) = 
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for any m,j. Thus e[ = Q^u^ ® u Xl ) = indeed. The other relations can also be seen 
by appropriately deforming the leftmost line from in the right diagram with the aid 
of the Yang-Baxter equation: 




Comparing the lines for pi in the left diagram here and the previous one, we find e2 + 
e' 2 + a' = e2 + a + d. Similarly the lines for p$ in the right diagram here and the previous 
one lead to es + e' 3 + b' + c' = es + b + c + d'. The proof is finished by noting d — d' = 
because of Qi(u m ® Uj) =0 for any m,j. □ 

As a corollary of Proposition 4.6 and (4.4), one has 

(4.15) Si(ut®p)=£i(p), 

which can also be verified by an argument similar to the above proof. 

Remark 4.7. Although the both and £, admit decompositions into the boundary 
and the bulk parts as in (4.1) and (4.13), these parts are not equal separately in general. 
Proposition 4.6 has also been proved by Mark Shimozono by using the technique known 
as katabolism (private communication). 

The energy £ n +i (4-12) and the row transfer matrix energy Ei (3.2) are related by 

Proposition 4.8. 

£n+l (p)-£ n+ i(T,(p)) = £7,( P ). 
For I = oo this coincides Lemma 4.5 with i = n + 1. 
Proof. We illustrate the proof for p = pi • • • <g>Pk with k = 3. Consider the diagrams: 



Ul Pi P2 P3 Woo Ul Pl p 2 P3 




Here the numbers above the vertices signify the n + 1 th non-winding number as in 
(4.14) with i = n + 1, and we have applied the Yang-Baxter relation to the line from 
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U[ . According to Lemma 4.4, £ n+ i(u l ® p) is equal to the sum of all the numbers in the 
left diagram. Similarly, £ n+ i(Ti(p)) is obtained from the right diagram as £ n+ i(Ti(p)) = 
d[ + d' 2 + d' 3 + a + b + c. The Yang-Baxter equation tells that + di = e\ + d[ for 
i = 1,2,3. Using these facts and (4.15), we obtain £ n+ i(p) — £ n+ i(Ti(p)) = £ n+ i(ui ® 
p) — £ n+ \{Ti{p)) = e[ + e' 2 + e' 3 , which coincides with Ei(p) in (3.2). □ 

4.4. Proof of Theorem 2.1. Theorem 2.1 is a simple corollary of (4.5) and 

Theorem 4.9. For any rigged configuration (^°\ (p^, r^), . . . , (f^ n \ r^™')) and the 
corresponding highest state p — p\ ® • • • ® pl under the KKR bijection, the associated 
ultradiscrete tau function (2.19) and the ultradiscrete CTM (4-1), (4-3) coincide. Namely 

(4.16) T i (p 1 ®---®p k ) = Pi(pi®---®Pk) (1 < i < n+1, 1 < k < L). 

Proof. Consider the embedding of p into V + (p {0) ) ® Bf L ' as p' = p ® \® L ' . The cor- 
responding rigged configuration is obtained from that of p by just changing p^ into 
^(o) y ^ j t j s eas jiy seen ^hat Tk i an( j p k i f or p' are same as those for p as long 
as 1 < k < L. Thus we understand them as associated with p' rather than p. 

Our proof is based on Propositions 5.1 and 6.1, which will be established in Sections 
5 and 6, respectively. Proposition 5.1 states that n satisfies the same bilinear equation 
(4.9) as pi. Combined with (4.8), it determines Tfc_i ; i, ■ ■ ■ , Tfc-Mi+i successively 
in this order from \i~k- i~k,ii ^~k.% |l<i<n + l}. Namely, the tau function on the 
NW corner in 

k 



Tk-l 


Tk 


Tk-1 





is fixed from those on the NE, SW and SE. Like p and p, the tau functions r and r 
are associated with p' and T oc (p'), respectively (see the beginning of Section 5), and the 
above diagram can be extended to a two-dimensional square lattice with the indicated 
coordinates. The square at (k,t) is associated with (^(^(p')))™^" . 

Consider the rectangular region on the lattice < t < t ,l < k < L + L' , where the 
tau functions for p' constitutes the top line t = of it. They are uniquely determined 
from the right boundary k = L + L' , i.e., {{tl+l 1 ,i(Tto{p')))™=\ I < i < t }, and the 
bottom boundary t = t , i.e., {(^(T^ (p')))™^ \ 1 < k < L + L'}. The coincidence 
of pi and n on these boundaries will be proved in Proposition 6.1 by taking to and L' 
sufficiently large. □ 

5. Bilinear relation for n 

Let TkA be the ultradiscrete tau function specified in Theorem 2.1 and (2.18)-(2.20). 
We define fk,d to be Tk,d with |sW| replaced by |sW| + [z^ 1 ^ in (2.20). In view of 
Proposition 3.5, this corresponds to the rigged configuration that has undergone the 
time evolution Too once. 

Proposition 5.1. The substitution pk,d — T~k,d an d~p~k.d = ^k.d solves the bilinear equa- 
tion (4-9). 

This section is devoted to the proof of Proposition 5.1 by a refinement of the approach 
in [26]. We invoke the free fermion construction of tau functions associated with g[(oo) 
[37]. For I e Z, set 

(5.1) a t (x) = (l\e H ^g\l), g = cxp(]T c ( ? ] {q^)) > 

(a,i) 
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where the notation is the same as eq.(2.3) in [37] except that 77 there is denoted by 
(7; here for distinction from (2.19). {pf^ here is not the vacancy number (2.7).) The 
operators il>{k) = X^ez V'jk 7 ; — Sjez _J ' are the free fermions. They obey the 

anti-commutation relations [tpi,ipj] + = [ip*,ipj] + = and [ipi,ipj]+ = Sij, hence tp(k) 2 = 
ip*(k) 2 — 0. \l) is the charge I vacuum of the Fock space. H{x) = J2 i>1 %i J2j e z ^j^j+i 
is the Hamiltonian with infinitely many time variables x = {x\, X2, ■ ■ •)• In (5-1), we 
associate each triple (c ^ , p , ) with the data {p^ , r| a ^ ) in the rigged configuration 
(A, (fj,^,r^), . . . , r^)). The sum extends over all the colors 1 < a < n and the 

rows 1 < i < %(°)). The tau function (5.1) is an ./V-soliton solution of the KP hierarchy 
with N = %«) + • • • + %(")). 

The time evolution of the free fermion is given by e H ^^p(k)e^ H ^ = e^ x ' k ^ip(k) and 
e H(xty*( k } e -H(x) = e -£(*,fc)^,*(fc) with ^( X;fc ) = Y, l > 1 x t k\ Consequently 

for x = eifr 1 ) := (/T 1 , ±/T 2 , i/T 3 , . . .). For z k := e^f 1 ) + • • • + e^ 1 ), the tau 
function is expanded as 

(5.2) (Ti(zk) = £ <ri(zk)v, 

(5.3, -^nw(C)'nf5^- 



(o,i) 



(5.4) A„ = 



n (v [a) -v {h) Ma {h) -a {a) ) 



where the sum (5.2) extends over the subsets C . . . , C independently. 
In (5.3), the product Il(a i) runs ovcr tne rows 01 the selected subset i^ a ) C In 
(5.4), ri(oi)<(6j) runs over tn e P a i rs °f sucn indices, whereas Il(a i) (6j) simply means 
the double product. A„ is the Cauchy determinant of the free fermion up to an overall 
power of and qf^ . It is derived by using the formulas: 



(m P W(q)\l)= £ p>q 



I 

q p 



m P i) • • • Hp m w{q m ) • • • r{q,)\D - n^-Pi)fe -go g 

iiij=i(Pi-«j) r=i 

Now we make a special choice of the parameters that further reflects the rigged con- 
figuration (A, {p {1 \r^),.. .,(//"), r("))). Fixing d E {2, . . . , n + 1} we set 

(a) (a) 

(5.5) = K^-5<i a) cM~ — ), q { ? ] = « (a+1) +^ (a) cxp(-^H, 



(5.6) cfV, 



(«)» 



5 ;W e xp(-^±^) if a e 
<5/ a ^ exp(~ M ' +ri ) otherwise, 



(5.7) ft =k w +5'jeM-—), 
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where 1 < a < n and e > 0. Here . . . , and 5\ a \ 8^ a \ 5" are e-independent 

generic (hence distinct) parameters such that 

(5.8) > • • • > k^" 1 ) > = > >■■■> 

(5.9) > 0, i (a) > 0, S'J > 0. 

Lemma 5.2. Set q = e _1 / e . In ffte limit q — > 0, £/ie summand oi{zk) v (5.3) of the tau 
function has the following behavior: 

o l {z k ) v = ^)+l^l+l^ ) l-'(k<- 1 >|-|^»l) (Xiy + (q)), Xv > 0, 

+ e{n^)) v = ^W-W^-W)^ + (q)), xl > 0, 

w/iere \v an d xl are independent of e. c(v, s) is defined by (2.20) with A = (Ai, . . . , A&). 

We denote by A a the relation a = lim e ^ +0 elog A under the ultradiscretization. 
It means that A = A q~ a + higher order terms in q for some leading coefficient A (j^ 0). 
(Aq still can depend on e as long as Aq ■— > although it is not needed in our case.) We 
let the relation A ~ B mean lim e ^ + o elog A = lim e _» + o elog£>. 

Proof. Let {{v[ a \ sf^)} be the subset of the rigged configuration r^)} corre- 

sponding to v = . . . , i/( n )) as in (2.20). We investigate the leading power of the 

constituent factors in (5.3). From (5.5)-(5.7) we find 

(i )J]c<«0gM ™ -^(l^l + I^D-l^l-l^l, 

(a,i) a=l 



% -iihm ™% (<, - i, i-i- w d. 



w n n - n n ^ - ^ 

(a,i) j = l Pj Pi i j = l Pi Pi 

(-) n (p^-p?)^ -?i a) )~nn^ a) -f5 o) ) 2 

(a,i)<(6,j) a=li<j 



UP 

a=l 



E(-min(!>\i/ a )) + 



(v) n (^ ) -^)- i -nn^ a) -^ l v i ^E min (- (Q) ,- (a - i) ), 

(a,i),(6,j) a=2 i,j a=2 

fa.fl K( ' "Pi • kW - Pi 



where = and the notation (2.3) is used. The contributions (i)-(v) sum up to 

-c{u,s) - - \v^\+l{\y^- r >\- \y^\). This verifies the leading power of a t (z k ) v in 
(5.10). Similarly, the one for oi(z k +£(k < - 1 ^ 1 ))^ is derived by including the contribution 
from (vi). 

The remaining task is to check the positivity and e-independence of the leading co- 
efficients \v and x'v We first illustrate them along o~i(zh)v In the right hand side of 
(5.3), we show the positivity individually for the constituent factors (i), (ii), (iii) and 
A„ =(iv)x(v) considered in the above. The leading coefficient from (i) is Y\, a ^ 8^ by 
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(5.6), which is positive due to (5.9). The leading coefficient from (ii) is 1 if I = 0. If 
I = — 1, it is given by 

Ad-i) (d+i) / ( a +i)\ e( ^ a) ^ 

K (d-1) 11 _ s (d) 11 ^ K (a) 



J a^d-l,d 



where the products on i and j extend over the selected rows in v (d ^ and v (d \ respec- 
tively. The symbol i{y^) denotes the length of as defined in (2.1). This is positive 
thanks to (5.8) and (5.9). The leading coefficient from (iii) with a fixed j is equal to the 
one from 

n K W - yr - ^ a+1 > 

1 1 Til \ , r(l) J-J- k-(I) - k(°) ' 

i o"q x J + o> 'q^i (<M) KW KW 



n>2 



It is positive by (5.8) and (5.9). The leading coefficients from (iv) and (v) are respectively 
equal to those in 

n n^v^ - *iv a) ) 2 n ^ {a) « (b) x« (b+i) - K ( a+i )))^ <a))£(,/b>) , 

a— 1 i<j l<a<6<n 
n 

n w-^ ] ^ a) - 4 a " i) ^°" l V 1 n ( K(a) - k^)-*^^^. 

a=2 i.i l<a,6<»> 

In view of (5.8) and (5.9), the coefficients are both positive apart from the same sign 
factor (-l)^i<a<i,<„ £( v a M v ) Thus the leading coefficient from the product (iv)x(v) 
is positive. 

For ai(zk + e(K^' -1 )), the leading positivity x'v > is proved similarly. The only 
necessary modification is to include the contribution from (vi): 

Y[ n a = M 1] - ^ a+1) Y {uia)) 



which is again positive due to (5.8) and (5.9). Finally, \v an d x'v are e-independent as 
they are rational functions of the parameters appearing in (5.8) and (5.9) only. □ 

Lemma 5.3. 

Co(^fc) > Tk,d, (T-l(Zk) ► Tk,d-1, 

<7o{z k +£(K y> )) >T M , CT_l(z fe +e(K V ; )) >Tk,d-l- 

Proof. For example we consider the UD limit 

(5.11) lim elo g( ro(^+e(« (lK1 )) = lim elogfVxU^+'^'Y 

where q — e _1 / e and (5.10) has been substituted. Lemma 5.2 furthermore tells that 
there is no cancellation in the c-sum here because of x'v > 0- Therefore the limit tends 
to max„{— c(u, s) — = Tk,d- See the definitions of Tfc(A) (2.19) and Tk,d (2-21). The 

other limits are confirmed similarly. □ 

Proof of Proposition 5.1. It is well known that ai satisfies the bilinear equation: 

(a" 1 - /?->,,(* + £(0 + eip- 1 ))*-!^ + e^ 1 )) 
+ (/3" 1 - 7 _ >o(* + eir 1 ) + e(7 _1 ))*-i(* + e^ 1 )) 
+ ( 7 _1 - a _1 to>(* + e(7 _1 ) + e{ a - 1 ))a- 1 {z + e^ 1 )) = 0. 
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This is derived by setting x = z + e^ 1 ) + e(/3 _1 ) + e(7 _1 ), x' = z and (I, I') = (0, -1) 
in eq.(2.4) M / in p956 of [37]. Setting 

a = « (1) , = Pk, 7 = oo, x = z k - 1 =£(/3r 1 ) + --- + £ (/3fe_ 1 1 ), 

we get 

(3 k a (z k -i + e(K (1)_1 ))cr-i(2;fe) 



(*k)<r-i(* fc -i + e^ 1 )" 1 )) + 5^e- A */Vo(«fc + e(^ (1) - 1 )) ( r_ 1 (z fe _ 1 ), 



where — k^ 1 ) has been evaluated by (5.7). In view of (5.8) and (5.9), the coefficients 
Pk, and 8' k ' here are all positive and e-independent. Moreover from Lemma 5.2, there 
is no cancellation of the leading terms coming from the two terms on the right hand side. 
Therefore by taking the UD limit lim e ^ +0 elog(-) of the two sides and applying Lemma 
5.3, we obtain 

(5.12) Tk-14 + Tk,d-i = max(Tfe j£ i + Tk-i,d-i, r k ,d + Tk-i,d-i - A&). 

This coincides with (4.9) with p replaced by r. Note that the range 2 < d < n + 1 for 
the both also match. This completes the proof of Proposition 5.1. □ 
Let us compare the results in this section with the similar ones in section IV of [26]. In 
[26], the tau function is supposed to fulfill the periodicity r k .d = Tk,d+n+i in the present 
notation. This led to a reduction condition (Prop. 4. 4 in [HHIKTT]) on each pair of the 
parameters (p\ a \ q\ ^ ) in (5.1), restricting the class of tau functions captured in the UD 
limit. In our approach, reduction conditions are bypassed by the special choice of the 
parameters (5.5)— (5.9) depending on the d that enters the bilinear equation (5.12) to 
prove. As it will turn out in Section 7.3, the ultradiscrete tau functions derived here 
cover all the solutions of the box-ball system. 

6. Asymptotic coincidence of n and p l 
6.1. Statement and its reduction. In this section we prove 

V 



Proposition 6.1. Given a highest path p with length L, set p' = p ® 1 ® • • • ® 1 and 

kg = L + L' . Then the equalities (1 < i < n + 1) 

(6-1) T k ,i{T^{p')) = p Ki {T^(p')) l<k<k , 

(6-2) rkUTlip')) = Pk ,i{Tl{p')) 0<t<t o 

hold if to ~S> 1 in (6.1), and if furthermore ko 3> Lto in (6.2). 

Combined with Proposition 5.1, it establishes Theorem 4.9 and thereby completes the 
proof of Theorem 2.1. Let {p^°\ {pS 1 \r^), . . . , (pS n \ r^)) be the rigged configuration 
corresponding to T^(p'). Without loss of generality we assume p^ < p$p < 
Moreover from the condition to ~> 1 and Proposition 3.5, we assume 

1 « rf» < r« < r« < • • • , 
(6 ' 3) r««rf i£tfUtf 

throughout this section. From Remark 3.4 and ko » Lto, the state T^(p') takes the 
form: 

a>l fc»l 



(6.4) {p') = u Xl <8> • • • <E> u Xl ® 1 <S> • • • <8> 1 <8>( ) <8> 1 <g> • ■ • <g> 1 



where A = (Ai, . . . , Xl) are the numbers such that p G B\ ± ® • • • ® B\ L and p e Bf is 
a highest path. 
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Lemma 6.2. Under the same condition as Proposition 6.1, the following relation holds: 

TkUTW)) ttwCZ^+V)) = Pk0 ATUp')) ~ PkoAT£\p')) (0<t<t o - 1). 

Proof. Suppose (p^°\ (p^, r^), . . . , ((jS n \ r^)) is the rigged configuration for Tj°(p'). 
From the definition (4.1) and the assumed situation (6.4), it is easily seen that Pk ,i(T^Ap')) — 
Pk ,i(T^ 1 {p')) is the number of balls with colors 2, . . . ,n + 1 contained in p, which is 
equal to \p^~>\. To calculate T ko , i ( T L(p')) ~ ^o.iC 7 ^ 1 ^)). wc a PP!y thc formula (2.23). 
Tk ,i(T^,(p')) is obtained by replacing A there with A U (l fe o-£) anc j r W pC 1 ) _ ^ Q _ 
by Proposition 3.5. Then the max contains fco only via min(A U (l fc(1_i ), v), hence 
one can let it be achieved at v = p^ by taking fco sufficiently large. Consequently 
T ko ,i(TL(p')) = min(AU(l fe o- i ), M ( 1 )) -min^ 1 ),^ 1 )) - \r™\ + (i - | +rf V (1) ) 

for any < t < t as long as fc > Lt Q . Therefore T koti (T^(p')) - ^(T^+V)) = \P {1) \ 
in agreement with Ph^M) ~ P/wC^+V))- □ 

By Lemma 6.2, (6.2) is attributed to t = to case. Thus the proof of Proposition 6.1 
reduces to showing (6.1), on which we shall concentrate from now on. 

Lemma 6.3. For 1 < k < L, p k ,i( T £(p')) = T kA T ^,(p')) = °- For L < k < k , the 
following relations hold: 

(6-5) PkAT^(p')) = Pk-L.M, 

(6-6) n,i(T£tf)) = Tk-L,i{p), 

where p £ Bf L is defined in (6.4-). 

Proof. For p kl i, the assertion is obvious from (6.4) and the definition (4.1). As for r k ,i, 
we use the expression (2.23) for Tj^(p') which corresponds to the rigged configuration 
(^,(^, r W),... I (^),rW)). 

(6.7) 7i(0= max{min(^i/)-miii(i/,i/)-|s|+T i (1) (i/)} £ C p (0) . 

According to (6.4), we have p^ = A U (1 L ). From Proposition 3.5, we know that 
r^ = p}pto + f- 1 ^, where fj 1 ^ is the rigging for p'. (This is also equal to thc rigging for 
p in Proposition 6.1 although this fact is not used below.) Thus to enters (6.7) only via 
\s\ = \is\to + |s|, where \s\ is to-independent . Fixing £ = (Ai,...,Afe) with 1 < k < L 
and taking t sufficiently large, we see that the maximum (6.7) forces the choice v = 0. 
This yields r feii (T^(p')) = ^ 1} (0) = for 1 < k < L, where the latter equality is due to 
Lemma 2.3. 

The maximum can be different from for L < k < fc , where we are allowed to take 
k so large up to fc depending on to- This corresponds to the situation (6.6), which will 
be considered in the sequel. To compute the right hand side of (6.6) by (6.7), we need 
to know the rigged configuration for p. In view of (6.4), it is obtained from thc one 
(/j(°',(/i( 1 ',rW),...,(/jW,f'W)) forT^(p') by replacing p^ with (1 L ') and thc rigging 
r^ with r^ = r^ — Y]^—i min(Aj, ^t^ 1 ^). Sec Lemma C.3. This amounts to changing 
\s\ in (6.7) to |s| — min(A, v) in the notation (2.3). Thus by setting £ = (l fe ~ L ), we get 



tu-lAp) = max {min((l fe L ), v) - min(^, v) - (\s\ - min(A, v)) + t\ X) (v)} 

k-L\ ,,\ _ •„/„ ,,\ _ U| , (1)/ 



= max {min(A U {r^ L ),v) - vam.{v,v) - \s\ + t> >{v)}. 

Since A U (l fc ~ L ) = (/4°\ . . . , p^), this is nothing but the expression of Tk^(T^(p')) by 
(6.7). ' ' □ 
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Thanks to Lemma 6.3, we may assume A = in (6.4) without loss of generality. 
To summarize so far, we have reduced Proposition 6.1 to (6.1) for p' such that p' G 
Bf k °. Resetting the meaning of p,p',p, L, L' and k 0l we restate it as 

Proposition 6.4. Let p 6 B® L be a highest path and ((l L ), (/i^, r^), (//"), r^)) 

be its rigged configuration with ^ < ^ < ■ ■ ■ ■ If L is sufficiently large and the 
condition (6.3) is satisfied, the equality 

(6-8) Tk,i(p) = Pk,iip) 

is valid for 1 < k < L, 2 < i < n + 1. 

A highest path p s Bf L satisfying the assumption of Proposition 6.4 will be called 
an asymptotic state. We have excluded i = 1 case since it is contained as the i = n + 1 
case of Tqo^p) which is also an asymptotic state. See (2.19), Proposition 3.5 and (4.2). 
The remainder of this section is devoted to the proof of Proposition 6.4. Our strategy is 
to express the both sides of (6.8) in terms of the quantities associated with the smaller 
algebra and invoke the induction with respect to n. Note that the induction allows 



us to use Theorem 2.2 with 1 < a < n — 1 and Theorem 4.9 for A 



(i) 

n-V 



6.2. Precise description of asymptotic states. The KKR bijection from rigged 
configurations to highest paths is known to be equivalent with the vertex operator con- 
struction [34, 35]. Here we utilize the notions in the latter formalism such as scattering 
data and normal ordering explained in Appendix D. In particular, we remark that a 
scattering data bi[di] (g>- • • (g> &at[c?at] 6 AS(B^) <g) • • -(g) Aff(B^) for an asymptotic state 
is normal ordered if and only if \i\ < ■ ■ ■ < ^n- 

Lemma 6.5. For an asymptotic state p, denote any successive tensor product components 
of the normal ordered scattering data by 

(6.9) 



B 



rfi 



Let the semistandard tableaux A and B be 



aia2 ...aij 



B 



bib 2 ...bit 



eBr 



2 < a\ < a 2 < ■ ■ ■ < ai A < n + 1, 
2 < h < b 2 < ■ ■ ■ < b lB < n + 1 . 



Then locally p has the form: 



(6.10) ■ • • 11 b tB ■ ■ ■ b 2 h 11 • • • 1 a u ■ ■ ■ o 2 aill 

Proof. Since p is an asymptotic state, we have Ib < I a- We divide the proof into two 
cases. 

Case 1. Assume Ib < I a- From the definition of the modes of scattering data (D.3), 
we have Is <C d\ — d 2 for asymptotic states. Therefore, the calculation of the vertex 
operator goes as (see around (E.2) for the explanation of <£>#) 

di —d 2 

$b(H ■■■lai A - ■■a 2 ai ■■■) 

di — d2—lB 

= h B - ■■ b 2 bi 11 • • • 1 Ti B (ai A ■ ■ ■ a 2 a\ • • • ) 

d\ — di 

= h B ■■ ■ b 2 b x 11 • • • 1 ai A • ■ ■ a 2 a x 

where T; is a time evolution of the box-ball system with capacity I career (3.1). 

Case 2. Next, consider the case Ib = I a = I- Let the energy function be H = H(B <g> A). 

Applying the definition of the mode (D.3) to (6.9), we have d 2 = l + rs + h and d\ = 
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I + 7^4 + h + H(B ® A), where r^, Tb are the riggings for A, B and h denotes the last 
term in (D.3) for d^ here. Since the asymptotic state satisfies the condition (6.3), we 
have tb < ta, leading to H < di — d 2 (=: A). If A > I, the proof is the same as Case 1. 
Therefore assume H < A < I in the following. Calculating the action of we arrive 
at the following situation: 



B- 



ai 



B ■ 



B' 



a, 



0,2 0,1 



where B' — 



l A &i& 2 



The diagram says that B 1 ® ai ~ a' ® (• • • ) under the 
combinatorial R. Let us show that a' = 6;_a- For the purpose, we first claim bi-\ < a;. 
In fact, suppose 6;_a > a-i on the contrary. We construct the pairs for B tg> A according 
to the graphical rule in Appendix B to compute H = H{B®A). We know that there are 
H winding pairs irrespective of the ways of making pairs. Since hi is weakly increasing 
with respect to i, we see that more than A + 1(> H) i's satisfy bi > a;. On the other 
hand, ai is the largest letter in A, therefore all the letters in B greater than ai have to 
constitute winding pairs, and we have seen that the number of these winding pairs is 
greater than H. This is a contradiction. Therefore we obtain bi-\ < a;. 

We have seen that 6;_a < ai 1 and we know that &z_a is the largest number in B' . When 
we construct the pairs for B' <S> a;, this fact means that 6;_a and ai form an unwinding 
pair. Therefore, the action of the combinatorial R is given by B' ® a; ~ a' <g> (■ • • ) with 
a' = 6;_a- By continuing the same argument, we arrive at (6.10). □ 

In what follows we use the notation explained in Section 2.1. 

Lemma 6.6. Suppose that Proposition 6.4 is true for A^}_i- For a rigged configuration 
((l L ),(Ai (1) ,r( 1 )),...,(^(™),r("))) with /iW = Oui,...,/ijv) and r^ = (n, . . . , r N ), as- 
sume fix < ■ ■ ■ < f-t j\[ . Then the corresponding scattering data bi[di] ® • • • <g) &jv[djv] £ 
Aff{B^ )®---®Aff(B%) is given by 

(6.11) &M = (X2, ■ ■ ■ , X n+ i), Xi = t\2j — T^_i ; j — 7"m^_1 + T M-l,j-U 

(6.12) d M = \fi [M ] I + r M + T^_ l n+1 - r$ n+1 . 

This lemma is shown without assuming that the scattering data b\[di] ® • • • <g> &Ar[<iAr] 
is normal ordered. 

Proof. From the arguments in Section 6.1, the assumption makes Theorem 2.1 for A^li 
valid. Then (6.11) is a corollary of Theorem 2.2 with a = 1. According to the definition 
(D.3), the mode c?m is given by 

(6.13) d M =VM+r M + H ( b J® b M +1) )- 

l<j<M 

On the other hand, combining (2.14) and (4.13) with i = n + 1, we have £ n+ i(bi ® • • • <£> 

& m) = Ei<j<m<M( min fe./ 1 '») _ H ( b j ® &™ +1) ))- (Since 6i <8> • • • (g) 6 M is A„_i-highest, 
the first term in (4.13) vanishes.) We know £ n+ i(bi ® • • • ®&m) = Pn+i(6i ® ■ ■ • ®6m) by 
Proposition 4.6. Moreover, since Proposition 6.4 for A^^ is assumed, we are allowed to 



use Theorem 4.9 to set ® ■ ■ ■ ® &m) = r Mn+i- Consequently, T M ^ n+1 is expressed 

as 

(6.14) rg n+1 = X! (min(^, Mro )-^(6 j ®6W+i) )) . 

l<j<m<M 

The formula (6.12) is a corollary of (6.13), (6.14) and the condition Mi < • • • < /ijv- O 

Given a rigged configuration ((1 L ), r^), . . . , (^ n \A n ^)) with fi^ = (/ii, . . . ,/ijv) 
and = (n, . . . , rjv), we introduce the numbers 

(6.15) k M ,i = min( ( u [M ] , /Z[ M ] ) - min( y u [M _ 1] , M[m-i] ) + r M + T~ M -i,i - t m^ 
for 1 < M < TV, 1 < i < n + 1. 

Lemma 6.7. Suppose that Proposition 6.4 is true for A^^. Let ((l L ),(fi^ ,r^), 
(fi^ n \ r^ )) be a rigged configuration for an asymptotic state. Set /j,^ — (fii, . . . , /i/v) 
with /ii < ■ ■ ■ < [An- Then the following relations are valid: 

(6.16) k M , n +i < k M , n < ■ < k M ,2 < k M ,i (1 < M < N), 

(6.17) k M ,i < k M +i,n+i (1 < M < N - 1), 

(6.18) fc M ,i - fcM,„+i = Mm (1 < M < N). 

Proof. By the assumption we may use Lemma 6.6. The scattering data bi[di] ® • • • <£> 
&Ar[c?Ar] considered there should be understood as a normal ordered one here because we 
deal with an asymptotic state and assume p,\ < ■ ■ ■ < /ijv. See the remark before Lemma 
6.5. From the definition (6.15), k M ,i-i - k M ,i = ^m-m-i _ T< M,i-i ~ T< M-i,i + T< M,i for 
2 < i < n + 1. This is equal to Xi in (6.11) hence nonnegative, proving (6.16). Summing 
this over 2 < i < n + 1 we get (6.18). Comparing (6.12) and (6.15), we have fcM,n+i = 
dm + \/j,[m-i]\- Therefore kM+i,n+i — ^m.i = ^M+i^n+i — ^M.n+i — Mm = djw+i — djvi- 
Since d^s are the modes of normal ordered scattering data, this is nonnegative, showing 
(6.17). □ 

Now we are ready to determine the precise form of asymptotic states from the associ- 
ated rigged configurations. 

Lemma 6.8. Suppose that Proposition 6.4 is true for A^li- For an asymptotic state p, 
let ((1 L ), (/U^> r ^)> * ■ ■ , (l^ n \ r^™))) be its rigged configuration and /xM — (/zi, . . . , ^ijy) 
with li\ < ■ ■■ < [In- Then p = p\ ® • • • ® pl £ Bf L is given by 



(6.19) Pk 



i k M ,i < k < k M ,i-i (2 < i < n + 1,1 < M < N), 

1 kM,l < k < fcM+l,n+l (0 < M < N), 



where fco.i = 0, fc./v+i,n+i = L. Namely p has the form: 

(6.20) 11 • • • ll(6i)ll H(6m)H ■ ■ • H(6m+i)H ll(M n • • • 11, 

where the segment (&m) S Bf MM (soliton) looks as 

kM,n+l ku,n kM,n-\ ' ' ' kus kuS-\ ■ ■ ■ ku,1 kjtf,l 

(6-21) I I | I | I I 

n+1, ■ • ■ ,n+l,n, ■ ■ ■ ,n, ,2,---,2 

Note that Lemma 6.7 guarantees that the regions of k appearing in (6.19) is the 
disjoint union decomposition of 1 < k < L. 
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Proof. By the assumption we may use Lemmas 6.6 and 6.7. In particular we use the no- 
tation Xi and du in Lemma 6.6. Lemma 6.5 tells that p indeed has the form (6.20). 
The segment (6m) has the left end at k' = du + |/"[m-i]| + 1 an d is arranged as 

n+l ■ ■ ■ n+1 • • • 2 • • • 2 with x% specified by (6.11). From the proof of Lemma 6.7, we 
find that k! = ku.n+i + 1 and x; t = ku.i-i — ku,i- Therefore it looks as (6.21). □ 

6.3. Evaluation of pi and Ti on asymptotic states. First we evaluate the tau 
function T ki of asymptotic states in terms of . 

Lemma 6.9. Suppose that Proposition 6.4 is true for A^^. If ((1 L ), (p^, f^), 
{p^ n \ r 1 ™))) is a rigged configuration for an asymptotic state with pS 1 ^ = (p,\, . . . , Pn), 
rW = (n, . . . , r^), (pi < ■ ■ ■ <pn), the associated tau function is given by 

(6.22) r k ,i = Mk - min(^[ M ], y u [M ]) - |r [M ]| + (k M ,i < k < k M +i,i), 
where < M < N,l < i < n + 1 and fc ,i = 0, /cat+i^ = L. 

Proof. From (2.23) we know 

(6.23) T kii — max {t(v)k — min(^, v) — \s\ + t\ X \v)}. 

Since s is the rigging attached to v and runs over the subset of A 1 " 1 that satisfies the 
asymptotic condition (6.3), the choice of v that attains the maximum must be of the 
form v = fi[ M ] for some < M < N. (We interpret /X[ ] = 0.) In terms of the notation 

(2.24) , we have t^(/X[ M ]) = t^\. In (6.23), the quantity in { } at v = P[m-i] an d 
v = /i[M] become equal if and only if 

(6.24) Mk - mm(p [M] ,p [M] ) - |r [M] | + rg. = (M - M - 1). 

This yields k = ku,i (6.15). Comparing the fc-dependence (M— l)k and Mk, we conclude 
that v = fi[ M ] gives a larger value than v = P\m-\\ if &M,i < k. Moreover we may use 
Lemma 6.7 by the assumption and therefore know that • • • < ku,i < ^M+i.i < ■ ■ ■ ■ Thus 
we conclude that the maximum in (6.23) is attained at v — p\u\ f° r k Md < k < k M+1 .i, 
where is equal to the left hand side of (6.24). □ 

Next we evaluate p k ,i for asymptotic states. 

Lemma 6.10. Under the same assumption as Lemma 6.9, pk,i for the asymptotic state 
is given by 

(6.25) pk,i=T k ,i (l<k<L,2<i<n + l), 
where the right hand side is specified by (6.22). 

Proof. By the assumption we may use Lemma 6.8, which specifies the concrete form 
of the asymptotic state as in (6.21). To evaluate Pk,i(p) (4-1), we count only the balls 
of colors 2,3, ... ,i in p itself and those of any color {2, . . . , n + 1} in the subsequent 
states T^ 1 (p). From Proposition 3.5 and (6.15), the positions ku,i in (6.21) changes as 
khi,i —> kM.i + Pm under the time evolution. Due to p\ < • • • < pn, there is no collision 
among the segments (solitons) (6m)'s in (6.20) under the time evolution. In view of these 
facts, the counting for pi-,i within the region kus < k < Am+m is done as 

M 

(6.26) Pk,i= ^2(k-k M >,i). 

M'=l 

From (6.15) and Lemma 2.3, this coincides with the right hand side of (6.22). □ 
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Example 6.11. The following figure helps to understand the counting (6.26). Consider 
an asymptotic state in which the M-th soliton is (6m) = 44332. Its time evolution takes 
the form: 



4 4 
T ' 

k = k M ,3 



3 3 2 



4 4 3 3 2 



4 4 3 3 2 



4 4 3 3 2 



t 

Here we have omitted ®, letters 1 and the other solitons for simplicity. Then the contri- 
bution to pk,3 from the M-th soliton comes from the balls within the frame, and their 
number is certainly equal to k — kM,3- 



Proof of Proposition 6.4- Due to Lemma 6.10 and induction on n, it now suffices to 
show n = 1 case of Proposition 6.4 to complete its proof. It is Lemma 6.6 that we 
started relying on the n — 1 case. But when n = 1, all the subsequent assertions are 
easily derived by only using Lemma 6.5 and the definitions of the scattering data and 
normal ordering in Appendix D. In particular, all the formulas are valid by setting 
t m\ = an d t m\ = — lM[M]l in agreement with the definition under (2.23). Thus (6.11) 
becomes bu = (^2) with x 2 = Hm, and (6.12) reads du — |M[M]| + r M- The definition 
(6.15) reads k M ,2 = k M ,i ~ Mm = min(/x [M ] , /z [M] ) - min(/x [M _ 1 ], H\ M -\]) + r M - Using 
the fact that tm < tm+i for normal ordered scattering data, one can directly verify the 
properties (6.16)-(6.21). By using them Lemma 6.9 is shown for n = 1, and (6.22) reads 
Tka = Tk,i + |M[M]| = Mk - min(/i[ M ],/i[ M ]) - |r[ M ]|. Finally (6.25) can be checked by 
substituting the above &m,2 into (6.26) with i = 2. This proves n = 1 case of Proposition 
6.4, therefore it is established for any n. □ 

Summary of proofs. We have finished proving Proposition 6.4. From the arguments 
in Section 6.1, it leads to Proposition 6.1. Combined with Proposition 5.1, Proposition 
6.1 proves Theorem 4.9 as explained in Section 4.4. Combined with (4.5), Theorem 4.9 
proves Theorem 2.1. 

In the course of these proofs, we have identified the three basic quantities by Proposi- 
tion 4.6 and Theorem 4.9. The tau function Ti (2.19) which is a piecewise linear function 
on the rigged configuration, the CTM for the box-ball system pi (4.1) and the energy £j 
(4.12) . We rephrase it as 

Theorem 6.12. For any rigged configuration (p(°\ (fi^\r^), . . . , (p^ n \r^)) and the 
corresponding highest path p\ ® ■ ■ ■ ® pl £ 'P+(t L ^)? the equality 

(6.27) Ti(pi <8> • • • ®pfc) = pi(pi <g> • • • ®p k ) = £i(pi ® • • • ®Pk) 

is valid for 1 < i < n + 1 and 1 < k < L. 

Note that the second equality (Proposition 4.6) has been shown even for non-highest 
states. The generalization of the first equality to them will be done in Theorem 7.4. 
Before closing the section we include a few immediate consequences. 

Corollary 6.13. For k = L, Theorem 6.12 becomes 

(6.28) n( Pl ®---®p L ) = pi{px ®---®p L ) = £i(pi ®---®Pl) = -c{n,r) - \p {i) \, 

where c(p,r) is the value of (2.20) at the "full choice" V(^ (a) ,s (a) ) = (^ (a) ,r (a) ), and we 
employ the convention \p^ n+1 ^ \ = as in (2.19). 
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Proof. For i = n + 1, the equality £ n +i(pi <8> • • • ® Pl) = — c (p, r ) is a consequence of 
the known relation between the charge of rigged configurations and the energy of paths 
[10, 22]. For i general, we find from (4.1) that p n+ \(pi ® • • • ®pl) — Pi{pi ® ■ ■ ■ ®Pl) is 
the number of balls with colors i + l,i + 2, . . . ,n+ 1 m pi ® • • • ®pz,. By the definition 
of the KKR bijection, it is equal to □ 

Remark 6.14. Corollary 6.13 tells that if A = the max (2.19) is attained at 

the full choice V(z», s^) = {^ a \r^). In particular, (2.23) leads to (//(")) = 
min(// a ),//° +1 )) - min(p( a+1 ),^ a+1 )) - |t> +1 ) | + r^ 1} (p (a+1) )- 

Now we are able to evaluate the conserved quantity E\ (3.2) for highest states in terms 
of the rigged configurations. 

Proposition 6.15. Letp E V + {p^) be the highest state corresponding to the rigged con- 
figuration (^ ) ) (/jW,rW),...,(/iW,rW)). Th en, its row transfer matrix energy Ei(p) 
(3.2) is given by Ei{p) = min(Z, £1^), which is E^ in (2.8). 

Proof. Combining Proposition 4.8 and Theorem 6.12, we have 

E t (p) = £ n+1 (p) - £n+i(Ti(p)) = t„+i(m (0) ) - < +1 (p (0) ). 

Here, by Proposition 3.5, T' n+1 (p,^) is obtained from T n +i(/i(°^) by replacing the rigging 

r^ with r^ = + min(Z, fi\^). This amounts to changing —\s\ in (2.23) (with 
a = 0, d = n + 1) into — \s\ — ^ min(/, On the other hand from Remark 6.14, we 

know that the max in (2.23) for A = /i^ is attained at v = /jp-h Therefore the difference 
7- n+ i(/i<°>) - < +1 (Ai< 0) ) is equal to £. min(/, ^. 1} ). □ 

Proposition 6.15 will be extended to non-highest states in Proposition 7.7. 

7. iV-SOLITON SOLUTIONS OF THE BOX-BALL SYSTEM 

As an application of Theorem 2.1, we present the solution of the initial value problem 
and iV-soliton solutions of the box-ball system. To cope with arbitrary states not neces- 
sarily highest, we first introduce in Section 7.1 an extension of the rigged configurations 
for such states, which we expect is equivalent to those studied in [23, 40]. We naturally 
extend the domain of the tau function to them. Generalizations of Theorems 2.1, 4.9 
and 6.12 to arbitrary (non-highest) states are presented in Section 7.2. Based on these 
results, we give the solution of the initial value problem in Section 7.3. In Section 7.4 
we derive several formulas for our tau functions in terms of the parameters that specify 
solitons. Together with (7.13), they yield the iV-soliton solution of the box-ball system. 
Our approach provides the general solution, which accommodates arbitrary number and 
kinds of solitons. A class of special solutions have been constructed earlier in [26]. 

7.1. n for non-highest states. For //W = (pf \. . . ,/i^ 0) ) £ (Z> X ) L , let p e B m <g> 
■ ■ ■ <g> B (0) be an arbitrary element not necessarily highest. Set 

(7.1) P = Pv ac ®P, 

(7.2) Pvac = (12 . . . n)® M " ® ■ ■ • ® (12) 8M2 ® l® Ml , 

where (12 . . . n) for example means 1 ® • ■ • <g> n G Bf n . The rigged configuration for p vac 
is given by 

(7.3) rc vac = ((l L °), ((1 L 0, (0 Ll )), • • • , ((1 L "), (0 L "))), 

n n 

(7.4) L a = ^(b - min(o, b))M b = ^ ~ a ) M *> (0 < a < «)• 

6=1 b=a+l 
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Thus L n = and ((1 L "), (0 L ")) actually means (0,0). The vacancy numbers pf (2.7) 
for the configuration ((f L °), (I Ll ), • • • , (1 L ")) of 

rc vac is calculated as 

n 

(7.5) 5 aA L -J2 C a,bL b = M a 

6=1 

for any j > 1. In (7.1), one can always make the state p highest by taking M\, . . . , M n 
sufficiently large. In fact, the choice 



(7.6) 



M a > m a+ i (1 < a < n) 



suffices, where m a denotes the total number of the letter a contained in the tableau 
representation of p. 

Let (jl, f) — (A , f )i • • • i (/* i f^ n >)) be the rigged configuration for the high- 

est state p. By the definition of the KKR bijcction, it "contains" rc vac (7.3) for p vac . By 
this we mean that (fi,r) can be depicted as follows (n = 3): 



/i 



(0) L 



1.(2) 








L 



Recall that (i^ is not limited to a partition, therefore it is not necessarily a Young 
diagram. Neither jj,^ has been depicted so. As mentioned after (2.5), any reordering of 
{{jl[ a \ f^)} for each a should be understood as the same rigged configuration. 

From the above rigged configuration (^°\ (p,^,f^), . . . , (^ n \f^)), we extract the 
data ( A i( 1 ),r( 1 )),...,( A i(™),r(")) by 

(7-7) M (a) = (a4 Q) )'=i U (1 L "), M (a) = (/4 a) )tr, 

(7.8) fW - (r> } + M a )<° =1 U (0 L «), r<«*> = (rf^l, 

for 1 < a < n, where l a — £(^ a ">). The shift M a in defining by (7.8) has been 
introduced on account of (7.5) and the algorithm for the KKR bijection, especially Lemma 
C.3. As the result, r^), . . ., (//W,rW) become independent of M 1 ,...,M n as 

they get large sufficiently. Therefore the data (fi, r) = (^°\ (fJ>^\ r^), . . . , (// n \ r^™^)) is 
determined unambiguously from p g 5 ( ) ®- • -<£>.B (o) by the prescription (7. 1)-(7. 8). We 

Mi Mi 

call (/i, r) the unrestricted rigged configuration for p, which we expect is equivalent to the 
one studied in [23, 40]. For highest states, it coincides with the rigged configuration under 
the KKR bijection, but in general (/Z^,// 1 ), . . . , /j,^) is not necessarily a configuration. 
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The vacancy number p) a > (2.7) can become negative. The rigging r^ G J} a is no 
longer limited to the range (2.6) but obeys the relaxed condition r- a ^ < p^ a) with 

some non-positive lower bound. We associate the tau function t^(A) (A C ^( a )) to an 
unrestricted rigged configuration (/i, r) by the same formula as (2.22). For A = /j,^) and 
p = pi <S> • • • <S> Pl, we will also use the notation r^A) = r kA = r^pi ® • • • ® p k ) as in 
(2.24). 

Example 7.1. Take n = 3 and consider the non- highest state p and the highest state p 
as 

p = 344 <gi 2 <g> 13 <g) 24 G B 3 <g) Bi ® B 2 O #2, 

P = Pvac ® P, 

p vac = 123123121 G £f 9 , 
where we have omitted (g> in p V ac- The rigged configuration (/t, r) for p is 













I— 1 




1 






























1 1 1 




1 













We have 

(Mi,M 2 ,M 3 ) = (1,1,2), (Lo,L u L2,L 3 ) = (9,5,2,0) 

according to (7.4). Thus the definitions (7.7) and (7.8) yield the unrestricted rigged 
configuration (/x, r) depicted as 



























Since p\ ' = —2, this is not a configuration. 
7.2. n — pi for non-highest states. 

Lemma 7.2. For any element p G B^(a) ® • •• ® ^ M <°); ^ Pvac,L a , (/2,f) and (^i,r) = 

(^'.(/iW.rW),...,^"'.^)) 6e as m (7.2)-(7.8). For a fixed A C the tau 

function associated with the rigged configuration (//, f) is decomposed as 

(7.9) n (\ U (1 L °)) = r 2 ( Pvac ) + L^(A) + Ti(A), 

(7.10) Ti(p vac ) = L Li - - ^ C a . b L a L b - L l 

l<a.b<n 

for sufficiently large Mi, M 2 , . . . , M n . Here Ti(p V3bC ) is the tau function for the rigged 
configuration rc vac (7.3). The last term in the right hand side of (7.9) is the tau function 
(2.19) associated with the unrestricted rigged configuration (p,r). 
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Proof. Let us write down the left hand side of (7.9) according to (2.19) and (2.20) as 
(7.H) 

t 4 (A U (1 L °)) = max|min(A U (1 L °), i> (1) ) - J V C a . b mm(u^,D^) - V |«(°) | - 

a, 6 a 

For Mi, M 2 , . . . , M„ sufficiently large, one has Lo 3> Li ^ ■ • • 3> L n -i 3> 1. In such a 
circumstance, one can show that the max can be limited to those C fS a ^ that contain 
(l La ) part entirely. Accordingly, we set 

= U (l 1 "), i»C M ( a ), 

|s(o)| = | s («)| + M a ^(z»), C t», 

taking (7.7) and (7.8) into account. Substituting these forms into (7. II) and using 
the formula (7.5) and min(t/ a ) U {l La ),v^ U {1 L »)) = L a L b + L a £(^) + L b t(v {a) ) + 
min(j»,^ b )), we obtain (7.9). The expression (7.10) is derived by means of (6.28). □ 

A decomposition parallel to (7.9) takes place also for p,. 

Lemma 7.3. Under the same setting as Lemma 7.2, set p = pi ® ■ ■ ■ ® Pl onrf take 
A = /itjjq in the notation (2.2), hence £(X) = k. Then for Mi, Mi, . . . ,M n sufficiently 
large, the following relation is valid: 

(7.12) p 4 (Pvac ®Pl ®---®Pk)= Pi(Pvac) + L±k + Pl ( Pl <g> • ■ • <g> Pfc). 

Proof. In view of p vac e Bf L ", the time evolution of p vac ®pi® ■ ■ - ®Pk under looks 
as follows (n = 3). 



L 

" M{ 



\L23123 ...12\ll .. 


. 11 




\ Pi(Pvac) 










L x k 



On the top row, the length Lq part is p V ac and the length k part is pi <E> ■ ■ ■ ® Pk ■ By the 
definition (4.1), pi(p vac ®Pi ® •• • ®Pk) is the number of balls with colors 2, ... ,i on the 
top row and all the balls in the SW quadrant beneath it. 

For Mi, . . . , M n sufficiently large, one has Lq ^> Mi ^> 1. Moreover from the time 
evolution rule in Proposition 3.2, the left segment within p vac with length Lq — Mi 
undergoes just a translation to the right by one lattice unit under T^. Thus this segment 
and the hatched region containing the balls arc entirely separated by the strip 11 ... 11 
of empty boxes with width Mi ^> 1. Therefore pi(p vac ® Pi ® • • • ® Pk) is decomposed 
into the contributions from p vac (trapezoid in the bottom left), pi ® ■ ■ ■ ® Pk (hatched 
region) and the parallelogram in the bottom. By the definition, the first two are equal 
to pi(p vac ) and pi(pi ® ■ ■ ■ ®Pk), respectively. The last one yields Lik because there are 
Li balls in total in the left segment in p vac with length L — Mi. □ 
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Now we give the generalization of Theorem 4.9 and Theorem 6.12 to arbitrary (non- 
highest) states. 

Theorem 7.4. For any state p = p\ ® ■ ■ ■ ® pl £ B (oj <S> ■ ■ ■ <S> B (o) , let (fi,r) = 

(//°), (/Z 1 ), r^ 1 )), . . . , (// n \r(™))) be the unrestricted rigged configuration, and let n be 
the associated tau function. Then the equality (6.27), namely, n(pi® ■ ■ ■ ®pk) = Pi(pi <8> 
• • • ®Pk) = £i(pi <8> ■ • ■ ®Pk) holds for 1 < k < L. 

Proof. The equality p, = Si has been already shown in Proposition 4.6 for any state, and 
we are only to show = pj. Since p vac ® pi ® • • • ®pk is a highest state associated with 
the rigged configuration (/^ U (l io ), (/iW.fW), . . . , f (">)), Theorem 4.9 tells that 

(7.12) is equal to (7.9) with A = Moreover it also tells that Ti(p vac ) = pi(p vac ). □ 

Combining Theorem 7.4 with (4.5), we obtain a generalization of Theorem 2.1 to 
arbitrary states. 

Corollary 7.5. For any element, p £ B^o)®- ■ ■®B^a), let (/i,r) = (fi (a \ (fj,^\r^), . . ., 
(^( n ), ?-( n ))) be the unrestricted rigged configuration. Then pk — ix\, . . . ,x n+ i) £ B < ) 
is expressed as 

Xd = TkA — T~k-l,d — Tk,d-l + T"fe-l,d-l 

m terms o/ ifte tow function t k.d = T d{{^ > ■ ■ ■ , M^)) associated with (p, r). 

7.3. iV-soliton solution. To simplify the notation we write A in place of in this 
subsection. We shall exclusively treat the states p = pi ® ■ • • <8>Pl £ B\ 1 ® • • • ® such 
that L is formally infinite and the boundary condition pk = u\ k is satisfied for k 3> 1. 
Under such a setting, the right hand side of the inequality (7.6) is still finite, therefore 
all the arguments in Sections 7.1 and 7.2 remain valid. 

Our solution of the initial value problem of the box-ball system is formulated as 

Theorem 7.6. For any initial state p — p\ <g> p 2 ® • • • £ B\ 1 ® B\ 2 ® • • • , Zet (/i, r) = 
(A, (z^ 1 ), r^ 1 )), . . . , r( n ))) be its unrestricted rigged configuration. Then the state af- 

ter the time evolution p'^p^®- ■ • = T^Ti 2 ■ ■ -Ti t {p) is expressed asp' k — (x\, . . . , x n+1 ) £ 
B\ k with 

(7.13) X d = Tk.d - Tk-l,d — Tk,d-1 + Tk-l,d-l- 

Here Tk.d = T~d{{\\, ■ ■ ■ , \kj) is the tau function (2.18)-(2.20) associated with (A, (/i*- 1 -*, r W), 

(//V (2) ), ...,(^«),r(™))) ; wherer'^ =r^ +^ t j=1 min(l j ,^). 

Proof. This is a consequence of Corollary 7.5 and Proposition 3.5. □ 

Let us evaluate the conserved quantity Ei (3.2) in terms of the data (A, (/z^, r^ 1 ), 
...,(//"), r("))). 

Proposition 7.7. For any state p = p\ ® pi ® • • ■ £ B\ x ® B\ 2 <g> • • • , let (fi, r) = 
(A, (fJ,^\ r^), . . . , {p^ n \ r^)) be its unrestricted rigged configuration. Then the row 
transfer matrix energy Ei(p) (3.2) is given by Ei(p) = J^j rnin(Z, fi^). 

When p is highest, this reduces to Proposition 6.15. 

Proof. Let p be the highest state (7.1) and let (fi,f) be the corresponding rigged config- 
uration. Proposition 6.15 tells that 

Ei(p) = £min(Z,/z«) =^min(«,( /U ( 1 )u(l^)),) = ii+E min ( Z " u ? ) )' 

3 3 3 

where we have substituted (7.7) into pf^\ On the other hand, due to Mi 3> 1 in (7.2) 
and the property (3.3), Ei(p) is decomposed as E[(p) = Ei(p vac ) + Ei(p). It is easy to 
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check E[(p vac ) = L x by counting the non- winding number using the graphical rule in 
Appendix B. □ 

Following [27, 26, 31], we call those states p of the box-ball system such that E[(p) = 
Y^j=i min(Z, i^p) /i-soliton states with amplitudes fi^\ . . . , p^\ Thus Proposition 7.7 
tells that any state of the box-ball system is an Zi-soliton state for some l\. More- 
over, Theorem 7.6 asserts that in the unrestricted rigged configuration (A, (p^\ r^), ■ ■ ■ , 
(^(™),r(™))), the A^l-t part (p^, (p( 2 \ r< 2 )), . . . , (//"), rW)) is the conserved quantity 
among which /i^ 1 ) provides the list of amplitudes of solitons. In the remainder of this 
section we set 

h=N, p {1) = (/*!,..., /ijv), = {r u ...,r N ), 
and rewrite the tau function in terms of the parameters that specify solitons. These 
parameters are equivalent to the conserved quantity (p, (p^ 2 \r^ ),..., (p^ n \r^)) as 
we will see shortly. The result yields the general 7V-soliton solution of the box-ball 
system, which supplements the special solution in [26]. 

From [30, 27, 31], it is known that A-soliton states in the A^ box-ball system are 
labelled with the affine crystal Aff(B| 2 ) ® ■■■ ® AS(B^ 2 ). The classical part 

-B,t 2 ® ■ ■ ■ ® -B,t 2 parametrizes the internal degrees of freedom of solitons. The affine 
part is incorporated in the integers ri, . . . , Vn, and specifies the positions of the solitons. 
Thus we start with any such data 

(7.14) h®---®b N eB} 2 ®---®B} 2 7 (n,...,rjv) e Z N , 

where we call each &j a soliton. Let (p, (p^ 2 \r^), . . . , (p^ n \r^)) be the unrestricted 
rigged configuration for bi ® ■ ■ ■ ® bjy. Without loss of generality we assume 

(7.15) H\ < ■ ■ ■ < jiN, i~i < rj if pi — pj and i < j. 

For any v C /i, let us express the A^}_ x tau function t^(v) associated with (p, 
(/i^ 2 \r^),. . ., (p( n \r( n ^)) in terms of bi, . . . , b^- We parametrize v C (p l7 . . . , p N ) as 
v = (pj 1 , . . . , pj M ) in terms of the subset J = {ji < • • • < j M } C {1,2,..., N}. From the 
array of N solitons bi ® ■ ■ ■ ® 6jv we extract an element in B^ 2 ® ■ ■ ■ ® B~ 2 by sending 
the corresponding components to the left by the combinatorial R as follows: 

B} 2 ® ® B} 2 ~ B} 2 ® ■ ■ ■ ® B~ 2 ® (• • • ) 

h® ®b N ~ bfj ® ■■■ ®b^ ] ® (•••)■ 

A caution is necessary about this notation. Consider for instance N = 3, M = 2 cases 
&i ® & 2 <8> b 3 ~ 6i ® 4 2) ® (0 for J = {1,3}, 
- & 2 1} ® ^ ® (0 for J = {2, 3}. 

(2) 

Obviously, the elements represented by the same symbol b 3 in the two lines arc not equal 
in general. In this way, b^ is uniquely determined only by further specifying J except 
a = 1. In what follows we will always take it for granted that J has been prescribed. 
From Theorem 7.4 for A^ 11 we know 

T P{v)=£ i {b^®..-®bW) (2< 4 <n + l). 
Applying the formula (4.13) for A^li to the right hand side we get 

^ ) M=E( & S+ & S + --- + ^)+^(^>--^C) (2<*<n + l), 

where tip = {b^\ 1 . . . , is the representation in terms of the number of tableau 

letters as in (2.11). The case i — 1 needs an independent derivation. We recall the 
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definition t^\v) = t^^v) — \ v\ given just before (2.22). Substituting the above formula 
with i = n + 1 to this, we find the result is unified into the single formula 

(7.17) rf V) = M-5>$h+- • • (! ^ ^ n+1 )' 

under the convention 



This is natural in view of the mod n structure of the indices in A n _ 1 . Similarly, the sum 
in (7.17) may well be written as b^J +1 + • • • + 6^ +2 . 

Now we are ready to express the tau function (2.23) associated with (A, (/U,r), 
( M <V (2) ),..., („(»>, rW)): 

(7.18) Tk.i — max{min(Aru , v) — min(z/, v) — \s\ + t- 1 * 1 (v)} (k > 1, 1 < i < n + 1) 

in terms of the solitons b\ ® • • • <g> 6jv and their positions n, . . . , r^. We parametrize ^ 
by </ = {ji, ■ • • ,3m} C {1, . . . , N} as before, and introduce the functions: 

(7.19) m ,iij) = min(A [fe] , N ) - rj - (6^ + ■■■ + 6« +1 + 6$) (j G J), 

(7.20) A 4 ( J) = 2 ^ min( W , Mm ) - # (&£> ® • • • ® 6^ 



l,m£J 
Z < m 



where min(Aru,/Uj) = X)m=i m in(A m , Mj) according to (2.3). (To simplify the formula, 
min(^,^ m ) has been kept as it is despite (7.15).) Substituting (7.17) into (7.18) and 
noting that mm{v,v) - \u\ = 2 1];,me,7,(<m mm (Mi> fj, m ) and \s\ = J2jej r j> we nnd tnat 
Tks is expressed as 

(7.21) r M = max { V - A,(J)| 

jc{i,...,jv} l-^ J 

for fc > 1, 1 < i < n+ 1. We introduce t^ = Tfc )Tl +i — |A[ fe ] | according to (2.18). Then by 
Theorem 7.6, the local states are specified by (7.13) and the time evolution T; is given 
by changing rj to rj + rnin(Z, /Zj), i.e., r)k,i{j) into - min(Z, fj,j). 

Using the formula (7.21), it is easy to evaluate the local state (7.13) explicitly for 
k S> 1 if Afe = 1 in this region and the condition (6.3) (without the super script "(1)" in 
the present notation) is satisfied. It yields the asymptotic state of the box-ball system 
well after the collisions of solitons. Omitting the derivation similar to Lemma 6.8, we 
give the final result: 

w 

(7.22) • • • 11 • • • ll(6i)ll 11(M ir^TT(6 M +i)ll 11(6jv)H • • • 11 • • • , 

where 1 g Bj and the symbol <g> has been suppressed. For each &m = (x 2 , . . . , x n+ \) 6 
B}1, (6 M ) G Bf^ M stands for the array 



n+1 . . . n+1 n . . . n 

In (7.22), the interval of adjacent solitons is given by w = tm+i — tm + S, where S is 
a constant independent of n, . . . , r^. Therefore if hm < fJ-M+i, we have w 3> 1 due to 
tm -C Tm+i- In case fj, M = fJ-M+i, we have 

(7.23) w = r M +i - r M + H(b M <g> 6 M +i) > H(b M <g> 6 M +i) 

because of (7.15). Here H(bM®bM+i) is the energy (2.14) for A^li crystals. It is known 
(cf. [27, 31]) that H(bM ® &m+i) is the minimum distance until which the solitons of the 
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same amplitude can get close. Therefore (7.23) is consistent with the fact that the tau 
function (7.21) constructed from the data (7.14) covers all the iV-soliton solutions. 

Our formula (7.21) possesses a structure analogous to the well known tau function of 
the KP hierarchy [37]. For each J, the sum J2jeJ Vk,i(j) is the superposition of individual 
solitons, whereas the quantity Aj( J) reflects a multi-body effect. A characteristic feature 
m f]k,i{j) (7-19) is that it contains in (7.16) rather than bj that appears in the 
asymptotic state (7.22). As for Aj(J), using the definition (4.11), it is "factorized" into 
the two-body function as 

(7.24) A l (J)= Yl Si(b<ff®b? a +1) ), 

l<(3<a<M 

(7.25) Si(b®c) = 2mm(l,m)-Qi(b®c) (b ® c e B} 2 ® B^ 2 ). 

Here ^f +1 ^ is determined by sending b^ in (7.16) to the left by the combinatorial R as 

6^ ® ■ ■ ■ ® ® 6 $f +1) ® ( ). 

Si in (7.25) is equal to min(Z, m) plus the i th winding number min(Z, m) — Qi{b®c). For 
i = n + 1, it has been identified as the two-body phase shift of the solitons labelled with 
b and c [27, 31]. Thus A, (J) can be regarded as a generalization of it to the multi-body 
phase shift for an arbitrary color i. 

7.4. Alternative forms of A^-soliton solution. We retain the notation in the previous 
subsection. The A-soliton solution (7.21) has been expressed in terms of the parameters 
in (7.14). Here we rewrite it further in terms of the scattering data (Appendices D and 
E): 

(7.26) &i[di] ® ■ ■ ■ ® b N [d N ] e Aff(^ 2 ) ® ■ ■ ■ ® AS(b} 2 ), 

(7.27) dj= rj + H{b k ®bf +1) ), & = 2'(Z>1). 

0<k<j 

Our task is essentially to switch from the position (rigging) rj to the mode dj. See (2.16) 
for the symbol 2 l . The mode dj here is a natural generalization of the one defined by 
(D.3). In fact, when b\ ® ■ ■ ■ ® bn is a highest element with respect to A n _\, one has 
b { p = 2»i and F(6 ®^ 1) ) = fj,j, hence (7.27) reduces to (D.3). The mode is transformed 
according to (A.l) under the combinatorial R. The afnnization of (7.16) reads 

Aff^ 2 )® ®AS(B} 2 ) ~ A$(B} 2 )®...®AS(B} 2 m )®(...) 

7 ' 28 h^} ® ® b N [d N ] ~ bV>[d$>] ®---® Ci4m } ] (• • • )• 

For the notation the same caution as for b^ is necessary as mentioned under (7.16). 
Applying the definition (7.27) to ^[d^] ® • ■ • ® ^'[rf-f] in the above, we find 

0<fJ<a 

where the notation is the same as (7.24) and we have employed the convention jo = 
and b^ = b a . The element b G Bf 2 in (7.27) is the A^_ x analogue of Uoo appearing in 

(4.12) for A [ n ] . By using (4.11), (4.12) and (2.14) for A^_ l crystals, this can be rewritten 
as 

dfj = r ja - £ n+1 (b^ ® ■ ■ ■ ® bfj ) + £ n+1 (b?> ® ■ ■ ■ ® bf~_l ] ) + E min (^= > ) - 

0<P<a 
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where min(^i 0) — Vjp- Taking the sum over a and using (4.13), we get 

M 

(7.29) E4: ) =E^+E 6 S-^ V + i(&g ) ®---®C)+ E nhn^M, 

a=i je,/ je,/ i</3<q<m 

where we have used + • ■ ■ + ^l+i = Mj- O n the other hand, from Corollary 7.5 we 
deduce 

M 

E^S + • • • + & = ^V) - - ^V) - ^i(^) + m a < * < n + 1), 

a=l 

where v = (fi^ , . . . , /j,j M ) as in the previous subsection. Since r} 1 ^ = £j for A^-i by 
Theorem 7.4, the right hand side here is evaluated by using (4.13), leading to 

(7.30) 

M 

E(S + 

= - E(^! + • ■ ■ + + w£> ® ■ ■ ■ ® O - ^pS? ® ■ ■ ■ ® C) + M- 

From (7.29) and (7.30), the quantity appearing in (7.21) is rewritten as 

M 

- Ai(J) = min(A [fc] ,i/) + E^ff + b & + • " + 

j£,7 a=l 

^•! )= 4?+^+ E min(^, M „)=r^+ £ S„ +1 (&^ ® &(f +1) ), 

l</3<a 0</3<a 

where S n+ \ is defined in (7.25). Thus we obtain 

M 

(7.31) r M - ^max {^(min(A [fe] ,/x j J-^ ) +6^ 2 + --- + 6H)} (1 < i < n+1), 

— ' "' a=l 

where the max extends over all the subsets J = {ji, . . . , Jm} C {1, . . . , TV}. Compared 
with (7.21), the expression (7.31) is formally free from the multi-body effect. It has been 
absorbed into the quantity <pf^ , which is a shifted mode. 

The formula (7.31) is most naturally presented in terms of the "principal picture" 
of affinc crystals rather than the conventional "homogeneous" one. To explain it, let us 
make a short digression on the principal picture in this paragraph. Recall that an element 
in the affinc A [ n ] crystal AS(B t ) is parametrized as (xi, . . . , x„+i)[(Z], where d G Z and 
Xi G Z>o are to satisfy x\ + ■ • • + x„+i = I. See (2.15). We naturally extend x, to 
i 6 Z by Xj+„+i = X{. Instead of (xi, . . . , x„+i)[<i], the element is also parametrized as 
Xi = 6i-i — 9i and d = 8 n in terms of an infinite sequence 9 = {Qi)i^% such that 

(7.32) ffi G Z, > ff u 6, = 9 t+n+1 + I for all i G Z. 

The correspondence between (x\, . . . ,x n+ i)[d] and is bijective. In fact, 6i = d — 
Xi — X2 — • • • — Xi for i > and #j = <i + xo + x_i + • • • + Xj+i for i < 0. We set 
Affp(Bi) = {8 = (6i)iez I (7-32)} and call the crystal structure induced on it the principal 
picture. Explicitly, it is given as follows: 

ej(6) = (6 t <5g +1) ), fj(6) = {9i + <5g +1) ) for 6 = (6 t ), 

where <5g +1 ^ = 1 if i = j mod n+1 and otherwise. If the right hand sides break the 
condition 6i-\ > 0i in (7.32), they are to be understood as 0. The combinatorial R is 
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especially simple in the principal picture: 

R : Aff p (B ; )® Aff p (B m ) — Aff p (B m ) ® Aff p (B,) 

(7.33) >— (^-50®(ei + Si). 

Here = Si+n+i = Si(8 (g) 9') is defined to be the color i two-body phase shift Si(b ® c) 
(7.25) for with b®> c e Bi ® B m , where 6 and c are specified by b = (#i_i - ) 
and c = (^_! - ^)r=i 1 - From (2.13), S l reads explicitly as 

(7.34) S i (6®6') = 2mm(l,m)-6 i + 6' i+n+1 - min {^ +fe -^ +fe _ 1 } 

l</c<n+l 

for 9 ®9' E Aff p (B / ) (g) Aff p (B m ). Observe the compatibility between (7.33) and (2.12). 
Actually for i = 0, the rule (7.33) on 6*o,#o disagrees with the changes of d, cf in (A.l) 
under the above mentioned identification 9q — d,9' — d! , which renders, however, no 
problem being merely the discrepancy in the normalizations of the energy function. By 
Affp we mean the crystal structure including the convention specified in (7.33). 9 is a 
generalized phase variable of solitons. 

Back to our iV-soliton solution, we restart with the principal picture of the scattering 
data (7.26): 

(7.35) 9 1 ® • • • ® 9 N G Aff p (^ 2 ) . . . g Aff p (&| 2 ). 
Accordingly, (7.28) reads 

Affp(B| i 2 ) g, . . . g, Aff p (5| 2 ) ~ AE P (B% ) ® • • • ® Aff p (^ 2 M ) ® (• • • ) 

(7.36) 

0i®---®0iv ^ ^ ) ®---®^ ) ® (•••), 

where, again, the notation is unambiguous only combined with J = . . . , j'm} as 
cautioned after (7.16). We set 6 ( "J = (^)iez and identify 9^\ with ^ in (7.31). 
OfJ e Aff p (B| 3 2 ) corresponds to (&JJ 2 , . . . , ^n+i)^] 6 Aff(B|£). Therefore we 
have efj ti = (f)fj - bf] a bf^ for 1 < i < n + 1. In this way (7.31) is simplified to 

M 

(7.37) r M = max {^(min(A [fe] ,MjJ-^)} (l<*<n+l), 

where the max extends over all the subsets J = . . . , j'm} Q {!)•••! -W} as in (7.31). 
Note that Oj"'^ = 0j" +1 + fj,j a is consistent with the time evolution rule in Proposition 
3.5 and T k ,i(p) = Tk^+^T^p)) indicated by (4.2). 

Finally we present an operator formalism that formally leads to (7.37) via the ultra- 
discretization. Let q be an indeterminate. Let A be the algebra over C[(7,<7 -1 ] generated 
by the symbols &(6), (6 G Aff p (_Bj- 2 )) that satisfy the commutation relations 

(9eAS p (B^),9'eAS p (B^)y. 

(7.38) *(0)**(0') = 

Here 6, 9' are related to 9, 9' by the combinatorial R (7.33), (7.34): 

(7.39) 9® 9' ~ 0' ® 0. 

(The commutation relation of ^ (9)^(9') and $>*(6)^*(6') arc not needed in the sequel.) 
We equip A with the time evolution T; (I G Z>i): 

IWJT," 1 = tf(T,(0)), T;**^- 1 = tf*(T,(0)), 

7}(0) = (fli + min(Z,m)) for = (9,) G Aff p (^ 2 ). 

T; is an automorphism of A since it commutes with the combinatorial R, i.e., T](0) ® 
T;(0') ~ Ti(9') ® T z (0) holds under (7.39). Obviously, T ( T m = T m Ti is valid. 
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For ieZ, let the bracket (•)$ : .4 — ► C{q,q 1 ] be the linear form on A characterized 
by the following properties: 

(7.41) (1), = 1, (X*(6))i = (X) u (V(6)X)i = q 6i {X) i for 9 = (fcjiez, 

where X denotes an arbitrary element in A. We shall write (Tj*XTf )i simply as (T l k X)i 

for any k £ Z. As an example, let 6 <8> <8> x 6 Aff p (.B| 2 ) ® Affp^- 2 ) ® Aff p (S^ 2 ). Then 
one has 

(T, + *•(*)) (*(0) + tf*(0)) (*( X ) + *•(*)) > 4 

= (T,*(0)tt(0tt( X )>. + <T,tt*(0)*(0*( X )>. + (T,tt(0)*'(0*( X )). + (T«*(0)*(fl**(x)). 
+ (T,**(0)**(0*( X )). + <T,**(0)*(0**( X )). + (T,*(0)*'(0tt"( X )>. 

+ (r i **(e)**(0)**( x )) r 

We need the following reordering of <8> <j> ® X by the combinatorial R: 

6®(t>®x^ ® 6' <8> x - ® X (2) ® 6»" ~ 6» O x (2) 0'^ X (1) ® ® </>'• 

See (7.36). As cautioned after (7.16), there are two elements x^ an d x^ that are 
relevant to x under the choices J = {2,3} and {1,3}, respectively. In terms of these 
elements, the above bracket is evaluated as 

I _|_ ^min(;,a)+ei _|_ min(i,6)+^ 1) _|_ min(i,c)+x^ 1 ' 

_|_ ^min(i,a)+min(;,fe)+e i +<^ i _|_ min(i,a)+min(i,c)+6l i +x- 2> _|_ ^min((,6)+min(i,c)+0< 1) + X f ' 
_|_ ^min(i,a)+min(i,6)+min(i,c)+^i+^»i+Xi 

From the commutation relation (7.38), the characterization of the bracket (7.41) and 
the definition (7.36), it follows that the tau function (7.37) associated to the scattering 
data 0i ® ■ • ■ (g> 0jv (7.35) comes out as the ultradiscretization: 
(7.42) 

k 

T k ,i = e limel g(n T^(*(6i) + • • • + (1 < i < n + 1), 

j=i 

where e is related to g by q = e _1 / c . The bracket is expanded into 2^ terms as in the 
above example (N = 3). In each of them, the list of the positions of "3/* specifies the 
subset J = {ji, ■ ■ ■ ,Jm} Q {1, ■ • ■ , N} for the relevant contribution in (7.37). The time 
evolution of the tau function r fej j(T;(p)) is obtained from (7.42) by further inserting the 
product rij=i ^X; 1 °f the automorphism (7.40). 

Unlike the tau function (5.1) for the KP hierarchy, A is not the Clifford algebra and 
it is not known to us whether the Laurent polynomial 
k 

<n T x'{noi) + • • • (*(M + ^*(o N ))). 

satisfies any sort of bilinear relations. However, the formula (7.42) is a most intrinsic 
way to present our ultradiscrete tau function. It synthesizes the principal features in 
the theories of solitons and crystal basis, i.e., the free-fermion like structure and the 
combinatorial R. 

8. Summary 

In this paper we have introduced the ultradiscrete tau function and exploited several 
properties related to the KKR bijection and the box-ball systems. 

In Section 2, Tj is introduced in (2.18)-(2.20) as a piecewise linear function on rigged 
configurations. The piecewise linear formula for the KKR bijection is stated in Theorem 
2.1. After a brief exposition on the box-ball system in Section 3, we have furthermore 
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introduced and £ j in Section 4. in (4.1) is the number of balls in the SW quadrant 
in the time evolution pattern of the box-ball system. £, defined by (4.12) and (4.11) is 
a sum of local energy function in the affine crystal. The fact pi — £i has been shown 
in Proposition 4.6. The two quantities provide analogues of the corner transfer matrix 
[1] in complementary viewpoints; pi from the box-ball system and £i from the crystal 
base theory. Theorem 2.1 is a consequence of the further identification r, = pi = £i in 
Theorem 6.12. Sections 5 and 6 are devoted to a proof of this fact. In Section 5, tj 
is shown to emerge as an ultradiscretization of the tau functions of the KP hierarchy 
(Lemma 5.3) and satisfy the Hirota type bilinear equation (Proposition 5.1). In Section 
6, n = pi is proved on the asymptotic states by induction on the rank (Proposition 6.1 
and its reduction in Proposition 6.4). These properties are enough to establish the claim 
T; = pi everywhere. Section 7 gives the generalization of Theorem 2.1 and Theorem 
6.12 to arbitrary (non-highest) states. As an application, the solution of the initial value 
problem in the box-ball system is given in Theorem 7.6. We have also included the 
formulas (7.21), (7.37) and (7.42) for general 7V-soliton solutions. Curiously, they are 
most elegantly presented in terms of affine crystals in the principal picture introduced in 
Section 7.4. 
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Appendix A. Crystals and combinatorial R 

The crystals B\ used in the main text are crystal bases of irreducible finite-dimensional 
representations of a quantum affine algebra U' q {o). Let us recall basic facts on them 
following [11, 12]. 

Let P be the weight lattice, {ai}o<i< n the simple roots, and {Ai}o<i< n the funda- 
mental weights of g. A crystal B is a finite set with weight decomposition B = U\ e pB\. 
The Kashiwara operators ej, fi (i = 0, 1, ■ • • , n) act on B as e% : B\ — ► B\ +ai U {0}, /j : 
B\ — ► B\— ai U {0}. In particular, these operators are nilpotent. By definition, we have 
fib = b' if and only if b = Bib'. For any b G B, set e,(6) = max{m > | e™& ^ 0} 
and (pi(b) = max{m > | f™b ^ 0}. Then we have the weight wt& of b by wtb = 
Et^(b)-e,(b))A t . 

For two crystals B and B', one can define the tensor product B ® B' = {b <g> b' \ b 6 
B, b' e B'}. The operators e"j, fi act on B <g> B' by 

e t b®b' if ipi(b) > e l (b') 
b®e t b' if ipi(b) < Si(b'), 



ei{b®b') = 
M ® ' I b®fib> \i<pi{b)<ei(b>). 



Here ® b' and b ® should be understood as 0. For crystals we are considering, there 
exists a unique isomorphism B ® B' B' ® B, i.e. a unique map which commutes with 
the action of Kashiwara operators. In particular, it preserves the weight. 

For a crystal B we define its affinization Aff(£>) = {b[d] d 6 Z, b e B} by ei(b[d]) = 
{eib)[d - Sio] and fi(b[d]) = {l<b)[d + S l0 }. {b[d] here corresponds to T- d af(b) in [12].) 
The crystal isomorphism B ® B' ^ B' ® B is lifted up to a map Aff(B) ® AS(B') 
AS(B') ® Aff(i?) called the combinatorial R. It has the following form: 

R : Att(B) <g> AS(B') — ► AS(B') ® AS(B) 
(A.l) b[d}®b'[d'} i — y ~b'[d'-H(b®b')]®b[d+H(b®b')}, 
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where b ® b' \— > b' ® b under the isomorphism B (g> B' B' <g> B. H(b ® 6') is called the 
energy function and determined up to an additive constant by 

C H(b®b') + 1 if 2 = 0, Mb) >e (V), Mb') >e (b), 
H(ii(b ® b')) = I H(b ® 6') - 1 if z = 0, ^,(6) < eo^), <M&') < e (&), 
[ H{b®b') otherwise. 

Proposition A.l (Yang-Baxter equation). The following equation holds on Aff(B) ® 
Aff{B') ® Aff(B"): 

(R <g> 1)(1 <g> <g> 1) = (1 <g> <g> 1)(1 (8) i?). 

We often write the map R simply by ~. The combinatorial R is naturally restricted 
to B® B'. 

In the main text we are concerned about the crystal B\ corresponding to the 1-fold 
symmetric tensor representation. We normalize the energy function so that 

max{iJ(6 ® c) | b ® c 6 B\ <g> B m } = min(Z, to). 

Under this convention one has min{7?(6 ® c) | 6 ® c £ Bj ® £> m } = 0. When I = m, 
the combinatorial R becomes the identity map on B\ ® but still acts non-trivially as 
i?(x[d] ® y[e}) = x[e - H(x ® y)] ® y[d + H(x ® j/)]. 

Appendix B. Graphical rule for combinatorial R 

Following [17], wc introduce a graphical rule to calculate the combinatorial R for An ^ 
and energy function given by (2.12) and (2.14). Given the two elements 

x = (xi,x 2 , ■ ■ ■ ,x n+ i) e B k , y= (2/1 , 2/2 , • • • ,2/n+i) e 
we draw the following diagram to represent the tensor product x ® y. 



Xl 




Vl 


X2 

























Combinatorial i? and the energy function for ® (with k > I) are calculated by 
the following rule. 

(1) Pick any dot, say in the right column and connect it with a dot in the 
left column by a line. The partner »' a is chosen from the dots which are in the 
lowest row among all dots whose positions are higher than that of » a . If there is 
no such dot, we return to the bottom and the partner »' a is chosen from the dots 
in the lowest row among all dots. In the latter case, we call such a pair or line 
"winding" . 

(2) Repeat the procedure (1) for the remaining unconnected dots (I — l)-times. 

(3) Action of the combinatorial R is obtained by moving all unpaired dots in the left 
column to the right horizontally. We do not touch the paired dots during this 
move. 

(4) The energy function H is given by the number of winding pairs. 

It is known that the results for the combinatorial R and the energy functions are not 
affected by the order of making pairs ([17] Propositions 3.15 & 3.17). For more properties, 
including that the above definition indeed satisfies the axiom, see [17]. 




By moving the unpaired dot (letter 2) in the left column to the right, we obtain 



1233 


<8> 


124 




133 


<8> 


1224 



Since we have one winding pair, the energy function is H I 1233 



<8> 



124 



1. 



For i £ Z n +i, the number of connecting lines that cross the horizontal level of the 
border between Xi and x i+ i is called the zth winding number. The energy function H is 
the n + 1 th winding number. The quantity min(7, k) — (i th winding number) is called 
the i th non-winding number. It is known that Qi(x <g> y) in (2.13) gives the i th non- 
winding number. By the definition, the winding numbers for x (S> y and y ® x are the 
same if x <£> y — y <%> x by the combinatorial R. 



Appendix C. KKR bijection 

In order to define the Kerov-Kirillov-Rcshctikhin (KKR) bijection, there are two dif- 
ferent ways. One is the original combinatorial algorithm [9, 10] explained here, and the 
other one is an algebraic version [34, 35] which will be treated in Appendix D. Although 
the both definitions are known to be equivalent, they work complcmcntarily in some 
aspects. In fact, we use the both definitions case by case in the main text. 

C.l. Definition. The KKR bijection provides one to one correspondence between the set 
of rigged configurations and that of highest paths. For a given rigged configuration 

(CI) RC=(( M f),(^ 1) ,rf ) ),---,(^ n) ,r5" ) )), 

we define the KKR procedure RC i — > p G B ( ) <S> ■ ■ ■ ® B ( ) <8> B ( > , which gives a 
highest path p. See Section 2.2 for definitions of rigged configurations, vacancy numbers 
and riggings. The data is called quantum space. 

Definition C.l. For a given RC, the image (or path) p of the KKR bijection is obtained 
by the following procedure. 

Step 1: For each row of the quantum space we assign the numbers from 1 to A 
arbitrarily, and reorder it as 

(C2) l* 0) = A 0) A 0) }- 

Take row fj,^ . 
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Step 2: We name each box of the row fj,^ as 



(C.3) 



Corresponding to the row ^ , let p\ be the array of l\ empty boxes: 
(C.4) 



a (0) 




(0) 


(0) 



Pi = 



Starting from the box a{ \ we recursively choose af' € ^ by the following Rule 1: 

Rule 1: Assume we have already chosen ^ € Let g^> be 

the set of all the rows of //W whose lengths w satisfy 

w > coZ(a^ -1 '), 

where the right hand side means the number of columns in /l/* -1 ) that 
are not located to the right of the box ot^ ^ . 

Let gi^ (c <?W) be the set of all the singular rows (<^$- rows whose 
corresponding vacancy number and rigging are equal) in the set If 
gi^ 7^ 0, then choose one of the shortest rows of gi % \ and denote its 



rightmost box by ol± \ If gf 1 — 0, then we take a\ 



(0 



a 



(n) 



0. 



Step 3: From RC, remove boxes a^°\ • • • , a^ 1 ^ chosen above, where j\ — 1 is 

the maximum k such that a^ fe ' ^ 0. After the removal, construct a new RC by 

Rule 2: Calculate the vacancy numbers p\ a) = £ 4 (a ~ 1} - 2E\ a) +E\ a+1) 
along the configuration after the removal. For those rows shortened by 
the removal, assign their vacancy numbers equal to the new riggings. 
For the other row, keep the original rigging before Step 3. 

Put letter j\ into the leftmost empty box of p\ as 



(C.5) 



Pi 



jl 









(°) 



a[°^ in this 



[iff in this order. 



Step 4: Repeat Step 2 and Step 3 for the rest of the boxes a 2 , a 
order. Put letters jk into empty boxes of pi from left to right. 
Step 5: Repeat Step 1 to Step 4 for the rest of the rows /j,^, M3°\ 

Then we obtain pk from /i^ , which we identify with the tableau representation of the 
element in B <o) . The image of the KKR bijection is given by p = pn ® ■ ■ ■ <S>f>2 <S>f>i- I 

The above procedure gives a map from rigged configurations to highest paths. Its 
inverse also admits a similar description. See Theorem 2 of [9]. 



C.2. Example of the KKR bijection. Let us illustrate a typical example of the KKR 
bijection. For a later convenience, we treat the single column type quantum space. The 
procedure for general quantum space is quite similar. 

Example C.2. We show that the following rigged configuration corresponds to a path 
p= 11112221322433. 



42 



,(°) 





(1) 
















X 


2 






3 





,(2) 







X 








oDo 



In the above diagram, we have specified the boxes to be removed by Step 3 with the 
symbol "x". Note that the boxes with "x" are the rightmost boxes of the shortest 
possible singular rows, and their column coordinates are increasing from the left to the 
right. We can remove three boxes at a time, thus resulting part of a path is 3 . Similarly 
we can proceed as 



(1 ) 
4 
4 



(l 12 ) 













X 


4 






3 



1 






X 
















X 


8 






3 








(i 11 ) o 



(1 ) 1 









X 






3 



E 

o □ 











X 


3 



oDo 



(I 9 ) 


2 










6 


X 


6 



Oxo 



oDo 



0X0 



(I 8 ) 







(I 7 ) 1 I x I x I x I l 
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(i 4 ) 



0' 





By removing all the boxes, we end up with 

p = 0®0®0®0®[2]®Q 

The following lemma is useful. 

Lemma C.3. Let p G V+(p^) and q G V + (i / ^) be the highest paths correspond- 
ing to the rigged configurations (fj,^,r^), . . . , {^ n \r^)) and (v^°\ (v^ , s^),. . ., 
{y{n)^ s (™))) ; respectively. Then the rigged configuration for the highest path p® q is given 
by 

(C.6) (m (0) U ^ (0) , U s'«), (//">, r* n )) U (i/M, a'W)). 

ifere (i/ a ),s ( a )) = {(f| a \ s/°^)} and i/ie rigging s ^ — (s/°^) is gwen 6y 



'(a) 



where pf^ 1 is the vacancy number (2.7) for (^°\ (fi^ ,r^), . . . , [pS n \r^)). 



Proof. Let q^ be the vacancy number for , (i/ 1 ) , s^ 1 ) ) , . 



s ("))). Then the va- 



cancy number p^ a ' for (C.6) reads p-^ = pj a ^ +9j°^ 
number — rigging) of the row (^ a ' ) , s/ "*) in (C.6) isp 



Therefore the co-rigging (:= vacancy 

'(<*) '(o) _ (a) 



-si -* with j 



(i/W, «("))). 



which is nothing but the co-rigging of the same row in (i/W, a* 1 )), 

Recall that the KKR procedure (Definition C.l) consults co-riggings to decide boxes to be 
removed from a rigged configuration. Therefore the above coincidence of the co-rigging 
means that the KKR procedure on (C.6) gives the path q when the part is firstly 
removed from U Moreover at this stage, the remaining rigged configuration is 
exactly (^ ),(M«,r( 1 )),...,(^"),r("))). □ 

Appendix D. Vertex operator formalism of the KKR bijection 

Here we give a crystal theoretic reformulation of the KKR bijection based on [34, 35]. 
The central notions are scattering data, normal ordering and the vertex operator. For 
illustrative examples, see Appendix E. 

D.l. Scattering data and normal ordering. We call elements of affine crystals &i[di]<8> 
• • • <8> b m [d m ] G Aff(_Bi 1 ) eg) • • • <g) Aff(Bi m ) scattering data. The number di is called the i-th 
mode. By using the combinatorial i?, scattering data can be reordered and the modes 
are changed accordingly. Given a scattering data s G AS(Bi 1 ) (g> • • • Aff(S; m ), define 
S m to be the set of such reordering as 

S m = {s' G □ ' AS(B ltr(1) ) ® • ■ • ® Aff(B, a(m) ) | s' ~ s}, 

where U' means the disjoint union over all the distinct permutations of (li, . . . , l rn ). For 
instance, if s = 234 ® 223 , we have 



I I 234 \ 7 



- ® 


223 


V 


234 




223 
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Note that in this case, the union over a is trivial as (I1J2) — (h,h) — (3,3), but S2 
contains two distinct elements since the combinatorial R is nontrivial as remarked in the 
end of Appendix A. 

For i = 2, . . . , m, let <Sj_i be the subset of <S, having the maximal i-th mode. Then 
we have 

(D.l) 0^5i CS 2 c ••• CS m . 



In the above example, we have S\ = { 234 ® 223 „}. We call the elements of Si normal 
ordered forms of s. In general the normal ordered form bi [di] ® • • • <g> b m [d m ] is not unique 
but the mode sequence d\, . . . , d m is unique by the definition and satisfies d\ < ■ ■ ■ < d m . 
Any element of Si is denoted by :s:. 

D.2. Maps CW,...,^"). Let (/zW, r^), . . . , (// n \ r^)) be an a£ ] rigged 

configuration. Pick the color a part (^ a \ r^). Here we simply write it as (/x, r). Namely 
/i = (fii, . . . ,Hm) is an array of positive integers and r = (rj), where rj is the rigging 
attached to the i-th row in /x of length zzj. For 1 < a < n, let Bi = Bf a+1 be the 

A^-a crystal in the sense explained around (2.17). Define the map C (a ^ among the A^l a 
crystals by 

C {a) : B lll <g> ■ ■ ■ ® B ilm — ► : Aff(B Ml ) ® • • • ® Aff(B„J : (1 < a < n) 
(D.2) 6i®---®6 m 1 ► :6i[di] <8> • • • ® M d m] : 

(D.3) <2 4 = r 4 + ,x,+ H(b k ®b\ k+1) ). 

l<k<i 

Here e £? Mi (j < i) is defined by bringing bj to the left by the combinatorial R as 

{bj ® ■■■® h-i) ® 6i ~ <g> (•••) 

under the isomorphism <8> ■ ■ ■ <S> S /Ji _ 1 ) <8> £? w ~ £? w ® (S^ ® • • • ® B w _,). Note that 
the choice (D.3) is compatible with (A.l). 

The map involves "A^ crystal" = {(«+ 1)*}. See (2.16) for the notation 

a 1 . The following suffices to define C^: 

(n + l) 1 ®(n + l) m ~ (n + l) m ® (n + 1)', Jf((n + 1)' ® (n + l) m ) = min(i, m). 

Since the normal ordering in (D.2) is not unique, is actually multi-valued in 
general. Here we mean by C {a \-) to pick any one of the normal ordered forms. is an 
operator that transforms elements of classical A„}_ a crystals to normal ordered scattering 
data by assigning the modes. 

D.3. Maps . . . , <|>( n ). Pick the color a and a — 1 parts of the configuration and 

denote them simply by = (/zi, . . . , fi m ) and ^ a ^^ = (Ai, . . . , Afc). Set £>; = Bf a+1 
and = -B ; . We define the map from the normal ordered scattering data in A n _ a 
affine crystals to classical A^]_ a+1 crystals: 

$W : : Aff(B Ml ) • • • ® Aff(B„ ra ) : -f ® • • • ® (1 < a < n) 

(D.4) 61 [dl] ® ■ ■ ■ ® 6m[dm] Cl ® • • • ® Cfe. 

From (D.3) and the fact that b\[di] ® • • • ® b m [d m } is normal ordered, we have < di < 
• • • < d m . Then the image ci ® ■ ■ ■ ® Ck is determined by the following relation under the 
isomorphism of A^l a+1 crystals: (We write T a d = \^f d € (B^ a )® d for short.) 

(T a dl O 61 ® T Q d2 ~ dl ® 6 2 ® • • • ® T^-^- 1 ® b m ) ® (a Xl ® a X2 ® ■ ■ ■ ® a Xk ) 

(D-5) 

~ (ci ® • • • <£> Cfc) ® tail, 
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Here we are regarding bi 6 B fli = B^ a+1 as an element of B' = B^ a by the natural 
embedding (2.17) as sets. The tail part has the same structure as (Xf 1 (g> 61 ® T^ 2 ~ dl ® 

• • -<H>&m) on the left hand side. In the actual use, it turns out to be (T dl (gxz^ 1 ®T d2 ~ dl ® 

• • • <g) a Mm ) containing the letter a only. (This fact will not be used.) 

To obtain c\ ® • • • ® Cfe using (D.5), one applies the ^l| l 1 2 a+1 combinatorial i? many 
times to carry (T dl ® 61 ® • • • <g) T^-^- 1 ® & m ) through (a Al ® a A2 ® • • • ® a Afc ) to the 
right. The procedure is depicted as 

a Xl a X2 a Xk 

a^ m 

a 

a 

1 1 

a 

a 

Cl C 2 Ck 

D.4. Vertex operator formalism. Define the A^p crystal element 
(D.6) p (n) = (n+ l)^ ® ■ ■ ■ ® (n + l)"'"' . 

Theorem D.l. TTie imagep of the rigged configuration (//°\ (/Z 1 ), r^), . . . , (// n \ r(™))) 
under the KKR bisection is given by 

(D.7) p = $«C (1) $ (2) C (2) • • • $ (n) C (n) (p w ). 

This is announced in [34] and proved in [35] . The theorem asserts that the right hand 
side is independent of the choices of the possibly non-unique normal ordered forms when 
applying the maps C^ 1 ', . . . , C^ n >. 

Set 

(D.8) = $(«+i)c(»+i)...$(")c(")(p(")) (0<a<n-l), 

which belongs to the A™ a crystal B-^ 1 ® • • • ® B-^ 1 . Thus p in (D.7) is 

Corollary D.2. for < a < n — 1, p^ coincides with the image of the truncated rigged 
configuration (p (a \ (// a+1 ), A a+1 ^), . . . , {^ n \A n ^)) under the KKR bisection. 

By the construction, C^(p^) is a normal ordered scattering data which is an element 
of an A„l a affine crystal. Then the map $( a ) produces an A„_ a+ i highest path by 

injecting the scattering data C^^p^) into the vacuum state a^ 1 ) ® • • • ® a^ 1 "- 1 . We 
call 3>(°) vertex operator in this sense. The construction (D.8) involves the family of 
scattering data and vertex operators for crystals of A^ C A^ C • • • C An \ It can be 
regarded as a crystal theoretical formulation of the nested Bethe ansatz due to Schultz 
[41]. 

Appendix E. Inverse scattering formalism of Box-ball system 

This appendix is an exposition of the inverse scattering formalism of the box-ball 
system mentioned in Section 3.2. We illustrate the calculations of scattering data, normal 
ordering and vertex operators explained in Appendix D along several examples. 



dm d n 



( a " 

1-1 

\ a - 

h 

< a 

di 

< a 
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E.l. Time evolution, scattering data and normal ordering. 



Example E.l. Consider the rigged configuration in Example C.2. We put many 1 £ B\ 
on the both sides of the corresponding path p = 11112221322433, and consider its time 
evolution under of the box-ball system. See Section 3.1 for the definition of T^. 



= 
= 1 
= 2 
= 3 
= 4 
= 5 
= 6 
= 7 
= 8 
= 9 



1111222211111133211143111111111111111111111111111111 
1111111122221111133211431111111111111111111111111111 
1111111111112222111133214311111111111111111111111111 
1111111111111111222211133243111111111111111111111111 
1111111111111111111122221132433111111111111111111111 
1111111111111111111111112221322433111111111111111111 
1111111111111111111111111112211322433211111111111111 
1111111111111111111111111111122111322143321111111111 
1111111111111111111111111111111221111322114332111111 
1111111111111111111111111111111112211111322111433211 



Here the length of the paths is 52, and t = 5 state contains the original path as l® 20 ®p<g> 
l® 18 . The following rigged configurations correspond to the above paths at each time. 



(I 52 ) 



38 
40 
43 



A 1 


(i) 


















7 






IE 


.+ 



(2) 



+ 4t 
f- 3i 
2t 















oDo 



The linear dependence of the rigging on t is in agreement with Proposition 3.5. The 
following is the list of all the normal ordered scattering data corresponding to each time 
t of the above paths. 



t = 0,l,2,3,4 



t = 5 



t = 6 



t = 7,8,9 



2222 4+4t ®L233] 11+3t 0L3^ 16+2t 



22 2 2 24 <g> 233 



2222 



24 



23 



20 



22 27 <Z> 2223 



2 (J 



22 



27 



223 



6« 


34 




334 




334 



20 



29 



2334 



2U 



29 



22 15+2t®l 223 lll+3t 



2334 



5+4t 



Compare this list with the above time evolution pattern. Each tensor product component 
of the scattering data corresponds to a soliton in the path. When the modes of the 
scattering data are well separated, the normal ordering is unique, and the corresponding 
path consists of well separated solitons that contain the tableau letters in the scattering 
data (in the reverse order), t ^ 5,6 are such cases. From the viewpoint of the scattering 
data, collisions of solitons happen when the modes get close and the normal ordering 
becomes non-unique, t = 5, 6 are such cases. See also Example 2.4 for the tau functions 
at t = 5, where Tks there is relevant to Tk+20,i here. 

Let us illustrate the derivation of the normal ordered scattering data at t = 5. At 
t = 5, riggings of /Z 1 ' attached to the rows of length 2, 3 and 4 are n = 23, r 2 = 22 and 
r 3 = 20, respectively. By Theorem D.l and (D.8), we know that p = ^^C^\p^), where 
C^\p^) is the normal ordered scattering data. It is constructed from the ^-highest 
path containing the letters 2, 3 and 4. According to Corollary D.2, p^ is the image 
of the KKR bijection of the following part of the original rigged configuration: 
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(1) 



(2) 



oDo 



Here /J- 1 ^ plays the role of the quantum space, and the KKR bijection is A 2 ' type with 
letters 2, 3 and 4. For example, if we can remove only a box from /i^ , then we have the 
letter 2 as a part of the path, whereas if boxes are removed from fj,^ , and /j,^ , the 
letter is 4. Removing the rows of //^ from the top, we obtain the A 2 highest path: 



(E.l) 



= h b 2 b 3 = 22 <g> 223 2334 . 



Assigning this with the modes according to (D.2) and (D.3), we get 



&i[di] b 2 [d 2 ] 63 Ms] = 22 223 2fi 2334 



25 ' 



To derive the mode d 3 = 25, for instance, we calculate X)i<fc<3 

H(b k ®b ( 3 +1) ) in (D.3) 

as 



22 ® 223 <g> 2334 ~ 22 ® 2223 ® 334 , 



where a ® 6 signifies the value of the energy function if (a ®b) = H. Since in (D.3), we 
have r 3 = 20 and /u 3 = 4, the mode is d 3 = 20 + 4 + + 1 = 25. 

To find the normal ordered scattering data C^(p^), we follow the procedure (D.l) 
and list the following sets: 



S3 = {22 



Si - 



25 



223 



222 



25 



23 



26 



23 3 4 25 ,L22J 25 



2334 



2222 



<S 2 = { 



22 



24 
25' 
- 24 

{ I 2222 



<8> 



23 



26 



8) 



334 



222 



2222 



2223 



25 



334 



26 



222 



2222 



23 



26 



334 



20 



2222 



25 " 
24 
25 
24 



2223 



2233 



233 







2233 



334 



26 ' 



34 



34 



26 



}, 



25 







34 



233 



26 



34 



26 



24 



23 



26 



33 4 26 ,L2222j 2 



233 



26 



341 }. 



The both elements in Si serve as the normal ordered scattering data in agreement with 
the previous list at t = 5. 

Example E.2. Here is a more intriguing example. 



(I 52 ) 



,(!) 



/( 2 ) 



38 
40 
43 



+ 4t 



5 + 3t 
10 + 2i 



oDo 



The normal ordered scattering data are listed below. 
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t = 0,1,2,3 
t = 4 



t = 5,6,7,8,9 



2222 
2222 
2222 



4+4t 



20 



233 



9+3t 



222 



222 



20 



20 



233 
2233 



23 



21 



22 



<8> 



2223 



<8> 


34 








334 


® 


34 


) 2334 




334 





34 

21 
21 
21 
21 
21 



13+2f 



22 20 ®L223J 21 ®|_2334J 21 



22 



12+2t 



223 



9+3t 



2334 



5+4* 



At i = 4, all 6 reorderings are simultaneously normal ordered. In a sense three solitons 
collide all together at t = 4. Compare this with the following time evolution pattern. 



1111222211113321143111111111111111111111111111111111 
1111111122221113321431111111111111111111111111111111 
1111111111112222113324311111111111111111111111111111 
1111111111111111222213243311111111111111111111111111 
1111111111111111111122132243321111111111111111111111 
1111111111111111111111221132214332111111111111111111 
1111111111111111111111112211132211433211111111111111 
1111111111111111111111111122111132211143321111111111 
1111111111111111111111111111221111132211114332111111 
1111111111111111111111111111112211111132211111433211 



E.2. Vertex operator construction of paths from scattering data. Here we illus- 
trate the action of the vertex operators <&W, . . . , introduced in Section D.3 (D.4). 
It is convenient to use the vertex type diagram to express the action of the combinatorial 
R. For example the following successive actions of the combinatorial R 



b' 



will be depicted by the diagram: 



Given a path p and an element b G Bi, one can carry b through p to the right by 
successively applying the combinatorial R as 



(E.2) 



b®p~p'®b', p,p'eB kl ®B, 



hi 



B 



few J 



under the isomorphism Bi <S> (B^ <8> ■ ■ ■ <S> B^ N ) ~ (B kl ® • • • ® B/, N )(£>£?;. As the result we 
get b' G Bi and another path p'. Actually, the only situation b' = ui (highest element of 
B{) will be encountered in our case, and the relation (E.2) will be denoted by &b(p) = p' ■ 
This is an elementary vertex operator. The previous ones . . . , defined by (D.5) 
are compositions of with several b. 



For example, to calculate <&r 



2334 



2334 



1233 



®5 



1123 



, the relevant diagram is 
1 1 1 



1112 



1111 



1111 



1 
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Therefore we obtain <E> | 2334 j (Jl"]^ ^ = 43321. Note that <!>(, has created one soliton 

labeled by the letters in b. 

In general, if bi [di] <g> ■ ■ ■ (£> 6 m [d m ] is a normal ordered scattering data, < J>( 1 ) defined 
by (D.5) is realized as the following composition of elementary vertex operators: 



(E.3) 



if 1 o$ bl oT? 2 ~ dl o 



■ot; 



d^-d„ 



1 o$. 



where f o g(p) = f(g(p))- The superscript "(1)" corresponds to that of /l/ 1 ). Note that 
for o = l, the effect of T a d in (D.5) is described by 7\ d = ($j-f|) d . 

In what follows we illustrate Theorem D.l and Corollary D.2. 

Example E.3. Take a path p = 11112221322433, which we have already considered in 
Example C.2 and Example E.l. From t = 5 case of Example E.l, the both sides of 



2222 



23 



334 



2222 



233 



34 



(E.4) 

serve as the normal ordered scattering data C^(p^). According to Theorem D.l and 
(D.8), the original path p is reconstructed as p = $( 1 )C( 1 Hp (1) )- This $W is realized, 
according to (E.3) and (E.4), as the following compositions of elementary vertex opera- 
tors: 



p = T? o $| 
= T i ° *[ 



2222 



ioT n 2 o$i ,o$ 



334 



ioT;o$, 



1 o <f> 

233 34 



(0"") 

(nf 4 ) 



It is easy to check p = 11112221322433 from these formulas. 

Let us illustrate Corollary D.2, which reflects the nested structure of the KKR bijec- 
tion. For <I>( a ) (D.5) with general a, the formula (E.3) is replaced by 

(E.5) $W = o $ 6l o ($ H )^ o ■ • • o o $ bm . 

Example E.4. We consider the same example as above. In the rigged configuration 
(see Example E.l), first look at the rightmost two diagrams which form an rigged 
configuration: 



/(2) 



oDo 



From , we set p^ = 4 according to (D.6) and obtain the scattering data (p^ ) = 
|~4~j , which is obviously normal ordered. From (E.5), the A\ highest path p^ = 
$(3)^(3)^(3)^ w i^h. letters 3 and 4 is constructed as 
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p (2) =$0O$H 

Taking the rigging attached to into account, we obtain the normal ordered scattering 
data 



C^(p^)=H 
Next we look at the following parts 

„.(!) 
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Then the A 2 highest path = $( 2 )c( 2) (p^) with letters 2, 3 and 4 is calculated along 
(E.5) as 



*B°*B o( *GD )3o<i> Gg4] 
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As a result, we have reproduced (E.l), which was the starting point of the previous 
Example E.3. Summarizing, the path p = 11112221322433 has been obtained as p = 

$(1)C(1)$(2)(7(2)$(3) C (3)^(3)^ 
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